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I. INTRODUCTION 


A. GENERAL CHARACTERISTICS OF INVENTORY MODELS 


Inventory models have been developed and used to deter- 
mine the optimal policy for running inventory Systems since 
ioilios The objective is usually t95 M@inimize the anntuai 
expected variable costs associated with running the inven- 
tory systen. Usually the costs are aggregated under three 


cost categories: 


1. The order costs depend on the freyuency of orders for 


2. The holding costs are those csostS associated with 


3. The shortage costs are those costs whick are incurred 


whenever the system runs out of stock. 


When a specific model 1S constructed and applied toa 
real world situation two additional goals tend to be sought. 
The first objective is to reflect within the model the rele- 
vant features of the specific type of real world situation 
to a satisfying degree. The second objective 1s to iimit 
the costs associated with applying the inventory model. 
These costs consist mainly of costs of data collection, Gata 
transfer,and computation. In general, a better result in 
attaining one goal implies a reduction in fulfilment of the 
other goal. 

Thais thesis will consider that special class of inven- 
tory models known as the singie period model. Such a nodel 
attempts to determine an optimal quantity of an item to buy 
and stock which will minimize the holding and shortage costs 


over the period's duration. Order costs are usually ignored 


because an order is assumed to be placed at the start of 


every given period. 


Be. DETERMINING THE SHORTAGE COST 


Of the three cost categories the shortage costs are 
usually the most difficult to express in jollar values 
because the data needed to evaluate such costs 1s seldon 
available ina readily useable form. Peeve deiy. 1S 
included as part of overhead costs or indirect costs. 

Assuming that shortage costs can be quantified, the 
usual approach is to associate the avent of 'no Gelivery' 
with a fixed number of dollars for each item not delivered 
immediately upon request. This can be refined by adding a 
time-dependent dollar cost which 15 a linear function of the 
length of the time perioi between the occurrence of the 
Shortage and the end of the period. An implicit assumption 
1s that the customer population is homogeneous. 

This thesis, considers the case where the customer popu- 
lation is not homogeneous. In this case sevaral different 
Major customers are assumed and there is a differert demand 
distribution and fixed unit shortage cost f9r each major 
customer. The unit shortage cost 1s assumed to be constant 


for eacn customer. 


C. THE NARF MODEL 


The specific model considered in this thesis is the 
Single period model presented by Y%cMasters [Ref. 1]. Phat 
molel is intended to provide supply support for overhauls of 
repairable aircraft systems done by Naval Airc Rework 
Facilities (NARF*s). The following assumptions apply: 

1. Within the modelled singie period, repair parts are 
needed for the overhaul of different end items. Each 


end item can be viewed as a different major customer. 


During the tepair of an era eit en, a demand fora 
given part occurs with probabitrey ee where p is 
greater than zero ,( and no demand will occur with 
probablatt yale) )re Demands are satisfied ona 
"first-come, first-served" basis. 

2. The repaic cycles for each type “or end item are 
nonoverlapping and one cycle iamediately follows its 
predecessor. There may be a different number of 
repair cycles per period for different end items. 

3. A fixed unit shortage cost is associated with a 
specific end iten. This unit shortage cost may be 
different for different end items. 

4. The guantity to buy for inventory is computed and 
placed in stock at the Leyinning of each quarter 
(period) and no replenishment during the quarter is 
allowed. It is assumed that the supporting Naval 
Supply Peenter Wisc) maintains the stock or will 
obtain it by the beginning of the quarter when it is 


" needed. 


Within this setting the expected shortage costs for the 


period are a function only of: 


1. The customer specific fixed unit shortage cost; 


2. The orijier in which demands occur during the period; 


3. The initial stock available to satisfy jemand. 


When a demand occurs, the amount demanded is one unit 
and the order in which each demand occurs depends on the 
probability distribution of the inter-demand times for each 
customer, where the inter-demand time is the time between 
successive occurrences of demand by that customer. The 
total demand by each customer during the period is binoni- 
ally distributed and is assumed to be stochastically inde- 


pendent of the demand of any other customer. 
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McMasters suggested an approximation formula for the 
expected shortage cost per period (two customers) (Ref. 1]. 
OO 
expected shortage ee eS Le iy YL £0x) Je 
a oe ae CX) CLA 
cost per period Wills y 
Om rOreadiscrete demarkis, 
y OO 
Tl, Wy Te fae 2 
ee F én a Be te FN } 
; (x= 7 ey 


i, 
where o- lag ee! 


Y = purchased quantity, initially available 
for delivery ; 

kK = X,+ Xi; 

X,= demand per period, customer 1 ; 

iS demand per period, customer 2 ; 

44> mean demand per period for customer 1 ; 

Az= mean demand per period for customer 2 ; 

2 ee iiesteolvone age COStmeror CUStONer 1 3 

alee ohORbadge ‘COst FOr eustomer 2; 
Pige= GenNSity function for x. 


P(x) = probability mass function for x 


The first term attempts to approximate the average unit 
Sieetage cOst. The second term evaluates the expected 
number of Shortages per period. The first term represents a 
convex combination of 7, and 1, weighted by jez and 4) « Its 
range of values is bounded by fz, anda %. A large value of 
Peeves the Valve Ofigthe first term towards #, whereas a 
large value of 4, drives it towards 7,. 

Several theses have expanded the work of McMasters. 
Slaybaugh [Ref. 2} analyzed the demand data and the rework 


schedule during 1980 for the TF-34 aircraft engine being 


reworked by NARF at the Naval Air Station Aiameda, 
Cal rrornid., Combining those data with the engine's bill of 
materials, he waS able to estimate the probability o£ 


replacement per repair action and conducted a parametric 
analyses to study the relationship between the shortage and 


SURPLUS COSTS. 
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Asselin [Ref. 3] appiied dynamic programaing win orager to 
expand the single period aodel to a muiti-~period model. ‘the 
development of a "good" formula Cor the expested Shortage 
costs per period inthe single period model is a basic 


building block for am Gificiene mull voce tego. 


D. PURPOSE 


The purpose of this thesis 1S to evaluate the approxina- 
tion for the shortage cost term suggested by 4ct#asters and 


to investigate alternatives to that approximation. 


E. PREVIEW OF ANALYSIS 


Chapter II shows that McMasters' approximation formula 
provides the exact expected shortage costs per period in 
case of an arbitrary number of mutually independent Poisson 
distributed customer demands. Because McMasters' approxima- 
tion formula does not provide the exact expested shortage 
costS per period in case of binomial distributed custoner 
demands, other approximation formulaS are presented in 
Chapter III. 

Simulation is then used in Chapter IV to VI to evaiuate 
the performance of the different approximation formulas. 
Chapter iV describes the simulation methodology. Chapter V 
contains Simulation results and discussion for the case or 
two customers each with binomially distributed demands. 
Chapter VI expands the analysis to more than two customers. 
The normal distribution is used to approximate tne sum of 
the binomial demand distributions and McMasters' approxina- 
tion formula is applied. This procedure is then analyzed by 
running Simulations. | 

A final summary with conciusions and recommendations is 
presented in Chapter VII. 


TZ 


II. POISSON DISTRIBUTED CUSTOMER DEMANDS 


Ae PRELIMINARY REMARKS 


Suppose that demands for each customer can be described 
by independent Poisson processes. In general, the Poisson 
distribution is often used to model infreguent randon 
events. For instance, the UICP inventory aodel of the 
US Navy for consumable items assumes the very low demand 
items to be Poisson distributed ( an average demani of 1 
unit or less per year). 

If the Poisson assumption is valid then the expected 
shortage costS per period are exactly equal to NcMasters' 
approximation and can be easily derived using conditional 
probability arguments. Because McMasters' approximation 
formula basically muitiplies an average of the unit Shortage 
costs by the expected number of shortages per period it will 
be named ‘monovariabie approximation formula’ (i.e. only the 


sum varies). 


B. COMPUTATION OF CONDITIONAL EXPECTED VALUE 


koss [Ref. 4 ppe-75-76] determines the conditional 
exmeected value §of X, given X;, + as. + Xp = X, wonere 
Xie ee ¢Xn ace independent Poisson random variables with 
Means My, 22+ «Hp +} 

First the conditional probability that X,= k given 


Le 


X, + oo + Xn = xX 1S ca leutated- 


1Ross computes the conditional expected value for the 
case of two Poisson distributions. His Poe of reasoniny is 
applied here to compute the conditiona tN value in 
Case, Or D2 Poisson distributions for any n 2 2. 


1s 


i 


. D : ex 
ikl: i) 


xX 
pt Aye tan -x{ 
Pie 


CX Se ee kf 


Pierre: z Ny ‘= x f 


> 
w 


: pax =kt PY X,+ ts een ae 
PLX 0 +X, 2x § 


Substituting the appropriate probability wmass functions 





yields: 

PL Xsk{ X71 +X, 2XS ay 
anh 22 (i, ee (eek) [| - _ 

- e’ tu Boe Cae edt Aw) Be nae fin) 
k! (ae eet Pog 





(x- k 
‘| uk Cae ae) . ) 
K 


x 
Cyt + Ap) 


cy 
> (i) aK (My po tiptlp ) 
¢ 


Lae Cli r + ta) Cyr cme 
2 (*) ( Ah ) ae ea aie 
k Lt tl “nt: + Mp 


Therefore the conditaohal “distran ition ee eco 


X, + e+ + X, = x is the binomial distribution with parane- 
J44, ; 

ters x and mer By symmetry it follows that the 

conditional distribution of X; given X,+ ... + X,=x 1s the 


z ; P : : : Mee 
binomial distribution with parameters x 21d ny 
£6) (1 =i) eres ie lle The associated conditional expectel 
values are ;: 


BC X71) y+ seein ; ea 


CmeveevaleOn, OF J bHE FORMULA FOR EXPECTED SHORTAGE COSTS 
PER PERIOD 


Given a shortage occurs, the demand generating the 
stockout is due to customer i with probability ae and the 
Eooaltie is 7, dollars. The expected penalty cost for the 
stockout is therefore 


Now the number of stockouts given x demands in the period is 


0 iP excy 
x-Y ast -x> Y ome 


The expected number of stockouts is therefore, 
OO 


ae 
OT Oa 1 Oe 
X= Y+1 


where P(X) denotes the probability mass function of X. 


Each of these stockouts results in an expected sost of 


atl 


Therefore tne total expected shortaye costs p2r period can 


be written as: 
oO 


i > 
= (Se (X= eee (x) 
Xx=Y+1 é 


The last expression is exactly the discrete analog to 
the approximation formula derived by UcMasters {Ref. 1] 
expanded to an arbitrary number of customer demand 


aretributrons. 
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D. SUMMARY AND CONCLUSIONS 


This chapter has shown that the monovarible apfproxima- 
tion formula provides a humerically exact value for the 
expected shortage cost when the aggregate demand process is 
the sum of arbitrary number of mutually independent Poisson 
processes. This formuia utilizes the closure theorem for 
the Poisson distribution, which says that tne sum of inde- 
pendent Poisson random variables is ajain a Poisson randon 
vaciable. 

Feller [Ref. 5 p.~.142] shows that the binomial distribu- 
tion with parameters (n,p) converges to the Poisson distri- 
bution with parameter “= ap as hw oo and p»0O in such a way 
that M=np remains bounded. So the Poisson distribution 
correspondS to a binomial distribution with large n and 
small p. 

Since the Poisson is a limit of binomials and since the 
formula for expected stockout cost is exact for the Poisson 
case, it seems reasonable that the formula might also be 
valid for the case in which demands from 2ach customer 
during a period are binomially distributed. Unfortunately, 
a single counterexample shows that the shortajye cost expres- 
Sion is not exact for every case in which the demands are 
binomial. To see this, consider the case in which there are 


two customers as follows: 


pee eo oa Probabilityeor stockout 
customer per perio demand per cycie cost 
1 1 Oo 2 100 
2 Z Or. 1 5) 1g,000 


We will show that expression (2.1) is noe Valid) fore nrc 
example. 

Let Xi, be the total demand in the period by customeri. 
Suppose X=X,+X,=2. Then 
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- 


PX + Xk, =2) P(X, +X,=2 | X,=k) 2 (X, =k) 


~ 
il 


0 
1 
= a P (X,=2-k) P(X, =k) 
nee = eeeCe =O) te Pik =1) P(x =1) 
2 1 2 1 


- one. 2 c 
ele Ca) (1 P.)P, 


0.1630 
ge, the Conditional probabilities are given by: 


P(X,=0, X,=2 
P(X =0 |X=2) = _P(Xy=0, X ee 


P (X=2) 
P(X, =0) P (X,=2) (1-p,) Pa 
ee a ee oie 
P (X=2) O26 30 
Pi eet (< 1) t) Pa (1~P2) P 
ee ait) _ fa) Pa (17 Pa) Pe | 0.9939 
1 Pa 2) 0. 1630 
@nd, the conditional expectation is: 
ECX, 1X=2] = OP (X,=0]X=2) + 1P(X,=1/X=2) = 0.9939 (2-2) 
Now let uS Compare (2.2) to the value given by expres- 
. eae oe, ‘ 2 --- fo. = 
ome (21) 3. 4, = Rp Pes aces Ee KC AZ |= eZ cannes 


1.6364. This is significantly larger than tha exact value 
Zee) = 

Thus we can see from the simple example that, for small 
ieee appropriate choices o£ the p;'s, expression (2.1) 
provides a poor approximation for E[X¢] X=x] when the total 
demand in a period from each customer is binomially distribd- 
uted. Consequently, the monovaciable expression for the 
total expected shortage cost will also be in error for this 
case. However, this example is extreme, and it mignt be the 
case that for some sets of values for the parameters the 


monovariable expression provides a reasonable approximation 


Wi 


for the true shortage cost. We evaluate this pecs rote 
a later chapter by simulating actual jietands Banaeeogtagce 
and comparing the monovariable expression to the average 
stockout costs computed over a large number of ceplications. 
Before we proceed to the development of alternative 
approximation formulas and to the SlMWtadrvomeg- seers eee 
first return to the simple example to understand another 
problem that results from the binomial demand situation. 
Suppose there is a single spare part to support both 
customers; i.e. Y=1. AS above, suppose that X=2. Now, what 
is the shortage cost? This cost is either $ 100 or $ 1,900 
depending on the seguence in which the demands occurred 
during the guarter. Below are three possibilities for which 
X=2 ( S corresponds to a demand during the period and * 


corresponds to no demand). 


S F 
Case 1 -=-= 422 -=-- ! ~aSe@.33 2. | =Seee——e oa | 
S £ S S 
Oona [cs=>5 | aos oe | 
S 
Case (23 “leel--..———— sa ] 
F S 
2 Nees oe Saas | 
in case 3, the stockout Cost is clearly $s 000-7 However, 


in cases 1 and 2 the stockout cost could be either $ 100 or 
$ 1,000 depending onwhen during a cepair cycle that the 
aqemands occurred. If, for example, it 1S assumed that 
demands occur at the end of a repair cycle, the stockeut 
cost in case 1 is $ 100, but in case 2 it would be either 
$100 or $ 1,000. If demands are assumed to occur at tne 
beginning of a repair cycle, the stockout cost ln case 2 
would be $ 1,000, but in case 1 it would be either $ 109 or 


$ 1,000 depending on what kind of rule is usei to determine 


18 


which customer getS anitem when two or more jemands occur 
Simuitaneously. 

Fe see, from this example, Lotro Ming.of the 
demands plays a key role in the determination of the 
shortage costs. We examine this issue Eurther in our Simu- 


fat nomecvaliation in the £ollowing chapters. 


ie 


III. DERIVATION OF OTHER APPROKIMATION FORAMULAS 


A. INTRODUCTION 


Before Simulations are introduced in the next chavter to 
evaluate the accuracy of the monovariable formula for the 
binomiai "real world" NARF Situation, two further approxina- 
tion formulas will be derived. Their accuracy will be eval- 


uated along with the monovariable formula. 


B. DERIVATION OF THE HYPERGEOMETRIC APPROXIMATION FORMULA 


Given two discrete mutually iniependent demand distribu- 
tions P(X) and P,(X,;) the probabrlity fom eacimepo-o. a 
demand vector <xX,,X,> can be evaluated from the product of 
their probability mass £unct lions; exw). 

Next in the conditional setting, i.e. Given a realiza- 
tion <x,,X,> with (X,+%,) > Y am interpretation Mormesameed 
drawings without replacements suggests an approximation 


based on the hypergeometric distribution. Let 


X humber of balis 1n urn (=x,+tx,), 


X,= number of balis with assoclated Unit Shottageueese | 


x, = number of balls with associa Codie weno t cade eeocsas 


1 *X,) = Vs 


S humber of shortages {= (x 


If s balls are drawn out of an urn containing x balls, 
x, Of which have associated unit shortage cost rm, and x, of 
which have associated unit shortage cost “, then, according 
to the formula for the expected valuS ot the ay pergeoncriic 
distribution, the expected number of x, balls in the sample 
is 
x 


——- ae a ee ee ee ee 


1 2 
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and the expected number of x, balls is 
es ; 


Multiplying each term by the associated unit shortage cost 


and summing them yields: 
x +x xX + xX ie (3.1) 


This expression represents an approxination for the expected 
Shortage costs per period given x,, x,, Y and the constraint 
that x, 


Because X, and X, are mutually independent random vari- 


ables, the unconditional approximation for the expected 
Shortage costS per period can be computed by multiplying 
each conditioned approximation with its associated prob- 


ability of occurrence and Summing these products, i.e.: 


Ty nN, 
2, Xt 


Ss . ; a | 
2 LD Amar (0,0) (So; % F oe o ed Cae” 
X20 X30 2 





This sum gives an approximation for the expected shortage 
cost per period. Unfortunately, it does not reduce to a 
Simple formula. The hypergeometric formula is programmed as 
the subroutine HYPGEO in the ;rogram provided as Appendix A. 
This formula corresponds to the first step of Mctasters' 


development of his approximation formula [Ref. 1]. 


C. DERIVATION OF THE CDF APPROXIMATION FORMULA 


In the conditional setting the hypergeometric approxina- 
tion formula does not reflect any information provided by 
the individual demand distributions. TO INcOLrpeLate this 
leeeidteron, sal daaitional weight can De introduced into the 


averaging operation in equation (3.1). 
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if (x,+x,) > Y , the guantity Y iS spilt neo seo 
subquantities Y, and Y, such that Xx,/x, = Y,/Y, . Suppose 
next that the density functions of f and Ge loon ave 
PLIOUL Cut. The areas represent the average height of the 
density functions in the indicated regions. 

In Figure 1a relatively snali fEract@on ope thom aon. 
ability mass over the region [0,x, ] fs locdted ein the waar 
(iy) Wiheted S fad relatively. large Eractlon more tie Sienmob- 
ability mass over the region [0,x,] is located in the region 
(Y, -X2J- 

The motivation for the CDF approxinatton Tommie la too 
the low fraction of probability Over (10x) Wie weed ec Sie 
fewer shortages for customer 1 and the high fraction over 
(Y5,X,] would result in more ShottageS  -oeeicustoucee 
Those fractions of probabidities “oven thom coon eee 
and (Y¥,,X,J ate incorporated as additional sveugnUSeeac ome ic 


hypergeometric approximation formula. 


B(X,) 
a 
— ne 
0 Y as 
B (X2) 
a | | 
Same ce oem (eae met ccm ome comes ces Tce meee eee fe ee ! oe 
0 Y X 
2 2 


Figure 1: Possible Density Functions for x and Xe 
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Let CDFt and CDF2 denote the cumulative distribution 
functions for X, and X, respectively. Then the additionai 


weights £, and £f, can be computed by 


CDF1(x,) - CDF1(Y,) 


a eee 
i CDF1 (x,) 
: CDF2(x,) - CDF2(Y2) 

a CDF2 (x,) ° 


These weights are then introduced as multipliers in the 
averaging operation or the hypergeometric approximation 


formula, i.e. eguation (3.1) becomes 


a Roe ee eee ns (3.3) 


The expansion to the unconditioned setting is identical 
to the procedure applied when deriving the hypergeometric 
formula at the end of the preceding section (1.e. Eyuation 
(3.2) ). The CDF approximation formula iS programmed as the 


Subroutine SUBCDF in the program provided as Appendix A. 


De. COMMENTS ON THE TWO DERIVED APPROXIMATION FORMULAS 


In case of the monovariable formula ther? iS ooniy one 
convex combination of the unit shortage costs computed, 
whereas by appiying the hypergeometric or CDF formulas a 
convex combination of the unit shortage costs for each 
possible demand vector must be conputed. Peet tgnt. bert ha t 
computing the convex combinations within the conditioned 
setting provides better approximations than the monovariable 
formula. 

In comparing these approximations under the Poisson 
distribution of Chapter II, the monovariable formula and the 
hypergeometric formula provide the same numerical results in 


the cases of mutually independent Poisson and binomialily 
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distributed demands when the latter is such that p,=p for 
aula. 

The suggested CDF formula attempts ~ftouretine tres. 
geometric approximation formula. However, it reguires the 
computation of the sukbquantities Y, and Y>. If Y, and Y, 
are not integers, a interpolation or the individual cunula- 
tive distribution functions has to be performed in order to 
compute the weights f£, and f,. An alternative approach would 
be to round Y, and Y, to the nearest Pnreqem ut opeinaees « 
both approaches cause distortions, therefore the CDF formula 
waS not evaluated by running simulations in the Poisson 
Situation. However, it is intuitively appealing and should 
be evaluated along with the monovariable and hypergeometric 
approximation formulas. The main question to be answered as 
far as the CDF approximation formula iS concerned, is 
whetker the incorporation of the distributions into the 


averaging operation is a step in the right direstion. 


Ee SUMMARY 


This chapter provided the arguments leading to the two 
further approximation formulas, called the hypergeometric 
approximation formula and the CDF approxination formula. 
Simulation is reguired to determine the accuracy to be 
expected from applying the three approximation formulas in 


the binomial "real world" situation of the NARF model. 


Zt 


IV. DESCRIPTION AND ANALYSIS OF SIMULATIONS 


A. PRELIMINARY REMARKS 


Ii the individual customer demands are binomiallv 
ZoeotEeputed with possibly different p;‘s and u;'s there is, 
based upon the knowledge of the author, no closed forns 
expression available which provides the exact value for the 
expected shortage costs per period. Therefore simulation 
was used in order to estimate the accuracy of the approxiaa- 
tion formulas.2 However, as with any Simulation, Bw ee 
provide only an estimate of the actual value of the expected 
shortage costs. The precision of the estimate gained from 
Simulation is almost exclusively determined by the number of 


Simulated periods (replications). 


Be PROGRAMMING LANGUAGE AND RANDOM NUMBER GENERATOR 


The simulation programs were written in FORTRAN WATFIV. 
WATFIV allows straightforward structured projram@Ming which 
improves the readability and verifiability of a progran. 
The WATFIV compiler provides a large set of error checks and 
gives, in case of errors, uncoded and precise error 
messages. A further reason is the availability of a WATFIV 
Specific debugging mode at the Naval Postgraduate School 


which eases testing of a progran. 


“In the early stages of the research done for this 
thesis the author had written Simulation programs for up to 
four customers whose demands, per period were assumed to be 
Poisson distributed. The Simulations suggested that’ the 
values of the approximation formula were in fact the exact 
expected Senta ge costs per period. Driven by this evidence 
the author finally succeeded in finding an anaiytical proof 
POE, tne Case Of a Poisson Setting. 
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The basic prereguisite for any Monte Caro sty pe 30: sami. 
lation is the availability of proper random naunzbers. The 
simulation programs use ‘LLRANDOMII', the Naval Postgraduate 
School Random Number Generator Package [Ref. 6]. aes 
package has been tested thoroughly and is well validated. 


It generates very good pseudo random numbers. 


C. SIMULATION OUTPUT ANALYSIS FOR MEAN ESTIMATE 


During Simulation many independent observations (repii- 
cations) of the shortage cost per period are obtained. 
Because, in the context of this thesis, the interest is 
limited to the mean of the Shortage costs per period, a 
straightforward output analysis can be applied.3 

The sample mean is used to estimate the mean of the 
Shortage cost per period. It is unbiased and has mininur 
Variance. Dividing the sample variance by the number of 
replications gives an unbiased estimate for the variance of 
the mean estimate. It represents the estimated precision of 


the mean estimate derived from Simulation. 


1. Notation 
C - scalar random variable, shortage cost per period 
C - sample average 
n - number of replications 
sd - standard deviation 


Var- variance 


A 
A hat indicates an eStimator, i.e. Var(C) is used 
to estimate Var(C). 


3The reasonin presented in this section is mainrlv 


extracted from meer 7}, as Eal” as Woutputee analy s=ae 
concerned. 
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' : S 
Var{C) = ----- ASG, = C)2 
ies 4 1=1 

Var(C) = --- Var (C) 
y, 
A 1 A ites 
sd(C) = ( -+ fax) | 
n 


By the Central Limit Theorem the distribution of the 
sample average of any random variable with a weli-defined 
mean and variance approaches more and more the normal 
distribution as the size of the sample becomes larger and 
larger. In general if there are more than twenty indepen- 
dent observations ina sample the normal distribution is 
assumed to provide a sufficient approximation. Therefore, 


when n is sufficiently large, the random variables 


e- EPee (Hae y at (C) 


> <a ce ce cee ee ee ee eee ee =e ePpeaRpaeamR SRP aD aa awe aman EP oe] aS as a 


sd (C) ; Var (C) 


ace distributed as a standard normal and a chi-~scuared 


distribution with (n-1) degrees of freedom, respectively. 


Both statistics are independent, therefore 


has a Student's t-distribution with (n-1) degrees of free- 
dom. Because the t-distribution converges to the normal 


Gistribution with increasing n it can be assumeji that 


is distributed as a standard normal random variable. 


Zy 


For the unknown true mean, a corfidence interval can be 


computed with ({1-x) as the preassigned contidence te ver: 


PIELC] wth Cz Le ed (C)} -1-« 


where Zax represents the standard normal deviate. 

To neease the precision of the estimated wmean of the 
shortage costs per period by a factor of two (1.e. reduce 
the estimated standard deviation of the mean estimate to 
half of its value) the number of replications has to be 
increased fourfold. In general, as far aS mean estimation 
through Simulation is concerned, the precision of the esti- 
Mate is directly proportional to the syuare root of the 
factor by which the number of replications gets increased. 
All simuiations run in connection with this thesis generated 


Sample sizes far beyond one thousand. 


De OUTPUT ANALYSIS FOR SIMULATED TWO BINOMIAL DEMAND 
SOURCES 


There are four different categories of ° parameters: 
humber of repair cycles per period, probability of one 
demand during one repair cycle, unit shortage cost, and the 
purchased quantity. Each of these parameters is a candidate 
for causing deviations of the numerical results of the three 
approximation formulas from the actual shortage costs per 
period. The tests of all three approximation formulas were 
restricted to two binomial demand sources in order to more 
eaSily identify those parameters which most aifect the accu- 
racy of the approximations. 

The output tableau of the simulation program and expia- 
Mations of its entries are provided as Appendix En The 
Simulation results presented in Chapter V for the different 
Settings display only those values from this tableau which 


are of intermediate interest. 


ae 


Ve. SIMULATION RESULTS FOR [WO DEMAND SOURCES 


A. PRELIMINARY REMARKS 


The program to perform the Simulations for two denzand 
sources iS provided as Appendix A to this thesis. The 
BOrlowinmg. SECtions of this chapter have the following 
Supucture: 

1. Purpose of the chosen parameter setting; 
2- Presentation of the Simulation results; 


3. Analysis of the simulation results. 


Be FIRST SIMULATION SETTING 


In the first Simulation setting the number of repair 
cycles and the initial seed are fixed for all runs. The p's 
and Y vary to analyze their influence on the precision of 
the approximation formuias. In Tables 1-4 the unit shortage 
costs are $ 100 and $ 1,000 for demand source 71 and demand 
source 2, respectively. In Table 5 the unit shortage cost 
for source 2 was reduced from $ 1,000 to $ 200 in order to 
analyze the effect of a lower variability among the unit 
shortage costs on the accuracy of the approximations. 

For each table the purchased quantity Y is a Specific 
multiple of the sum of expected demands per period fron 
both demand sources. The value of Y is always rounded to 
the nearest integer. However, Y was increased by at least 
one unit in the sequence of Tables 1-4 in ordjer to examine 
tiemenernence OL 04a Ghange in Y om the accuracy o£ the 
approximations. 

The simulation results presented in this and _ the 
following section are arranged in tableaus consisting of 
blocks of six numbers which refer to specific parameter 


values. 


29 


Shortage- *% 
cost SD/MeEn Nn -2 y 7 


__DEVMON © ___DEVHYP_6____DEVCDF_®| 
The meanings and units of the abbreviations are the same 
as described in Chapter IV, section )D. In the tableaus, P1 
and P2 denote the p values for the binomial distributions 
foc demand source 1 and demand source 2, respectively. The 
maximal sum of both demands per period is 9 units for 
Tables 1-5. 


1. Parameter Setting for Tables 1-4: 7” 

Number of time units 20 

Number of repair cycles(1) 4 

Number of repair cycles(2) 5 

Unit shortage cost(1) $ 100.- 

Unit shortage cost (2) > 1000.- 

Initial seed for all runs: 8 150449 

*The observed sample mean for shortage costs per 


peciod derived from the Simulation. 


. ,>The observed standard deviation for the sample nean 
divided by the sample mean provides a measure (in percent) 
of the precision of the simulation. 


6Deviations of the values of the three cifferent 
aioe ma ton formulas from the simulatior result in 
percent. 


values, £or see 


“In order to avoid having very large 
=0.1, Y=0 was chosen 


Tables 2-4 for the case of P1=0.1 and P2 
ror that case in Table 1. 


“For all simulations in this sasect lene the same uiiaer ae 
seed was used in order _ to reduce the variability orf _ the 
Simulation results. Differences in the performance of the 
approximation formulas between the tables are therefore 
mostly due to different levels of Y. 
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because Y 
Therefore no deviations could be computed. 
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The following summary contains the average absolute 
deviation (in percent) of the approximation formulas from 


the simulation results: 


level monovar. hypergeon. Cur 

Or Y rormula formula formula 
Table 1: 0.85 mean 7.4 5.3 1320 
Table 2: hean 1. 6 11.4 igs 
Table 3: 1.15 mean Zesac2 HS 5 205 
Table 4; 1.30 mean S05 24.7 30a.5 
peo he) : 1.30 mean 14.8 ere 19.0 


with reduced unit shor- 
mlge COSt Variability 


ine this setting the hypergeometric formula shows 
Slightly more accurate results than the monovariable 
formula. The CDF formula shows significantly less accurate 
results than the two other formulas. 

Tables 1-4 reveai that for all three approximations the 
accuracy of the approximations decreases as Y increases. 

The accuracies of the approximations acre highly depen- 
dent on the ratio of the two unit shortage costs of the two 
demand sources. If both unit shortage costs were the sane, 
the problem is reduced to approximating the expected number 
of shortages per period. The higher the variability of the 
unit shortage costs, the less accurate is the approximation 
of the expected shortage costs per period. Lomo Lustrate 
this property, the variability of the unit shortage costs 
was set lower in Table 5 than in Tables 1-4. As can be seen 
when Table 5 is compared to Table 4, the change of unit 
shortage cost(2) from $ 1,000 to $ 200 reduces the inaccura- 


cies of the approximations to haif of their former size. 


3 


Although Table 1 represents simulation cresultsS with a 
high variability of the two unit snortage costs, the averaye 
absolute deviation for the sixteen runs 1S 7.4 percent with 
a Maximal absolute deviation of 16 percent for the monovari- 
able approximation. For the hypergeometric formula the 
average absolute deviation is 5.3 percent with a maximun 
absolute deviation of 13 percent. These results suggest 
that the monovariable and the hypergeometric formulas 
provide an acceptable accuracy when Y is less than the wean 
of the sum of both demands. Theretore, further analyses 
will be restricted to settings where ¥ is eéjyual or greater 


than the mean of the sum of demands. 


C. SECOND SIMULATION SETTING 


In this setting the number of repair cycles per demand 
source were equal. To analyze the effect of smaller 
increases of Y the three levels are now: = mean, Y= 1.19 
mean (i.ee. only 10 percent increase), and Y= 1.25 mean. The 
values of P1 and P2 were taken to be 0.1, 0.3, 0.5 and 0.3. 


Number of time units 10 
Number of repair cycles(1) 10 
Number of repair cycles (2) 10 
Unit shortage cost(1) > 000g 
Unit shortage cost(2) S$ 100.- 


Initial seed for all runs: 40150 
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5. Analysis of thewsornularron  esuteoe 


—_— = << eee ee oe ee eee Se ee ee es ee ee ee ee ee ee es ee ee es ee ee ee eee ee ee 


The average absclute deviations (in percent) of the 


approximation formulas from the simulation results are: 


level monovar. hypergeon. CDE 

or y formula formula forauzula 
Table 6: Y = mean 4.9 a6 16320 
Table 7: Y = 1.10 mean vS5 Sone 17a 
Table 8: Y = 1.25 mean 11.17 4.4 ZO 2 


Tables 6-8 indicate again the relatively strong depen- 
dence of the accuracy of the approximations on the size of Y 
relative to the mean. Although the variability of the unit 
shortage costs is hiyh, the monovariable and the hypergeore- 
tric formulas seem to provide acceptable accuracy even for 
the case where Y = 1.25 mean. Within this very special 
Setting all three approximations indicate tnat the accuracy 
decreases Slightly as the difference between P1 and P2 
increases. However, Tables 1-8 do not indicate in all cases 
that the precision of the approximations increase as the p 
values decrease, aS might be expected from the fact that the 
Poisson distribution is the limiting distribution for the 
binomial distribution when n becomes large and p_ becomes 


small. 


D. THIRD SIMULATION SETTING 


This setting is a Variation on the Qeerecerscecng. In 
order to get a more general survey about the performances of 
the approximation formulas the initial seeds anld all paranme- 
ters but the unit shortage costs are different for each run. 
The unit shortage costs are restricted to $ 100 and ¢ TC GgG. 


On2 row refers to one Simulation run. 


J 


1. Table 9: Runs with Kiaely Varying Parameters 
N1 N2 P1 P2 P11 Pi2 iMerorlD DEVNON DEVHYP DEVCDF te 
Seeomeee ace O20, 1000 100 2.6 3° 259 es 230 506 
Se 40.2 0.7 100 1000 3.4 4 999 0 = 12 -1 
peeoeeoeo O24 100 1000 4.2 5 255 14 24 ae 
eee. 250.6 «€©6©6100 1000 55.0 8 100 236 136 31a 
fae Oso 0.1 1000 100° 2.5 3 777 a2 10 35 
6 5 *0.4 0.9 100 1000 6.9 7 654 30 14 ag 
eee Ons 0.1 1000 100 95.8 7 379 ae ae = 
Peo. oo.s 1000 100 3.4 5 111 84 52 111 
ceeplbeeG.7 0.5 100 71000 4.7 5 7 98 118 i> 1 
ewe oruso 1000 100 5.8 6 825 245 216 325 
Meeemeoeo Oy 100 1000 2.9 3 900 a2 20 36 
pee Oe.2 02-9 1000 100 5.0 5 693 -24 3 =.) 
peeoeso 0.5 10008" 100 5.7 7 387 17 4 33 
Seer os | O27 100°>1000 8.0 9 951 —1 -4 =] 
cum OReOTG 1000 100 6.6 8 291 ee 34 15 
Paden runs 20,000 replications 

2. Analysis of the Simulation Results: 

The average absolute deviations (in percent) are; 
DEVMON = 74.7, DEVAS — 6) =, 1, and DevepDF = 117.5. These 
error deviations are relatively large. However, a reduction 


§ 100 and 


absolute deviations by about 


Pomeeveom Vartavility of the unit shortage cost to 
$ 200 would reduce the average 


one-kal£& demonstrated in Tables 4 and 5. 


l2Abbreviations in Table 9: 


N1 = humber of repair eee. 5 

N2 =] hum oer OL Bonet cycles (2 

Pi - p parameter 

P2 - DPD parameter (2 

Pl1i - unit shortage oo 

Pi2 - unit shortage cost (2 

MU - expected sum of demand per period 

x —- purchased quantity. 

SSObadmcrthitial Seed Wsed for result in Current row. : 
DEVMON, DEVHYP, and DEVCDF represent the deviations of 


the approximations from the simulation result in percent. 


Sai) 


E. FOURTH SEMULATIONSS er ViNG 


This simulation setting 1S Similar to both the second 
and st ned sSetban Ge As in the second setting, the number of 
repair cycles 1s the same for both demand sources. AS in 
the third setting, the parameters have row by cow comparable 
values. A comparison between corresponding rows in Table 9 
and Table 10 should indicate whether the performances of the 
approximation formulas improve if the lower number of repair 
cycles is increased such that both demand sources have the 
same number of repair cycles. The rcatio Y/4U is about the 


same aS in Table 9. 


1. Table 10: Identical Parameters _ as Tabie 9; But N1=N2 


Bile eee ee Pil Pi2 MU Y SEED DEVMON DEVHYP DEVCDF 13 
3 3 0.8 026 100071005 ale 12 1 34 
4H & 0.2 0.7 100 10005532055 wea 11 73 14 
5 5 0.6 0-45 5 100 100GNS 20 Cyne =a) > = 
G9 =~ 022 026 0d oe 9 12 Ass 13 at Zo 
4 7 0.6 0.1 1000 1005) Soo ial 7 = 4 
6 6 O.4 0.9 1001000 57223 c - waho sis 31 
7 7 0.8 0.1) 100055 100 Sic23> 7 ware 6 me) 6 
4 4&4 Q0.6 0.4 1000 100 4.0 5 123 g 0 17 
6 96) 0.7 0.5 00 1000 =2 ca. = + “2.0 
? 7.0.8 02-6 1000 100 »926>5 70 5350 3 a2 21 
3 3 0-5 0.7) 9) 100" 1000R ou eco 7 = 17 
4 7 -Q.2 029 1000 100075 7 Soe =17 15 —4 1 
9 9 0.6 0.3 1000 100 }o. lige, 1 =) v 
171.11 0.10.7 100 1000 8.8 10 449g cl ail 5 
7 


7 0.4 0.6 1000 10005 750 ames -1 7 = 14 
Fach runs 20,000 replications . 


13The abbreviations are the Same as for Table 9. 
(see preceding footnote) 
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Peat ysis Of (thes samulation Results: 


The average absolute deviations (in percent) are: 
DEVON = 7.5, DEVHYP = 4.2, and UEVGRE =— 10. 1. The siluwu- 
dation results establish that the monovariable and the 
hypergeometric approximation formulas provide better accu- 
racy when the number of repair cycles is the same for both 


demand sources. 


Fe FIFTH SIMULATION SETTING 


The increase of the lower number of repair cycles of 
Table 9 such that both demand sources had an equal number of 
repair cycles resulted ina dramatic improvement in the 
accuracy of all three approximation formulas. This improve- 
ment could be caused by two facts. The first 1s that the 
approximations improve when there is an identical number of 
repair cycles, as established in the second and fourth sinmu- 
iation settings. The second is that an increase of the 
parameter nis one condition for the convergence of the 
binonial lise see sLreyel to the Poisson iass tL 1 Dutaome. 
Therefore the improvement could be caused simply py the 
increased number of repair cycles. 

Table 11 displays results of Simulation runs with the 
humber of repair cycles of one demand source of Table 10 
further increased. The ratio Y/MU is about the same as in 
Table 10. Comparison of Tables 9,10, and. 11 provides an 
indication of both the impact of increasing the n values and 


having them equal in value. 


3) 9) 


1. Table.11: Continuation of Tabic 0 eevee ore 


increased 


Niet! eee Pil Pi2 MU ¥Y SEED DEVHON DEVHYP DEVCD 
3 5 0.8 0.6 1000 ~T0098 525 cue (Zz 43 W 31 W 77 W 
5 & O.2 0.7 100, 1000 ese 2 —2 | isa = Sree 
5 7 0.6 0.4 100 1000 9535e5 7) S232) -ioa =] Wo = 23a 
9 10 0.2 0.8 100 1000 9.8 13 9822 =119) 27 = 2k 
9 7: 0.6 0.1) 1000™ 1009 oes 17 rag = Ze 
6 7 0.4 0.9 100 10000 S27 9 Gee =] 1 Sia ae 
8 7 028 0.17 1000> 100 Wei coe 0.1 “= S tOsaa eee 
G 6 0.6 0.4 1000 “100042 3c me 26 W Ie 41 W 
6 8 0.7 0.5 100 1000 8.2 5° 29950-1530 — it oe 
7 9 0.8 0.6 1000 S100 i ieaice 17 W 11 =W ein 
5 3 °0.5 027 1000000” So= se 20 W a BN 
7 8 0.2 0.9 1000 100 Seno > yi 8 I 51 NS 2S 2 
119 0.60.37 1000" ViItG Sai = es a aS 1 i 
22 11 O21 O27 “100 > 100G 3-3 I ieene 15 W ay 22 
9 7. 0.4 026 1000” 100 2 S3o Waeey - ion <7 5 =e 


Each run: 20,000 replications 


2. Analysis of the Simulation Results: 


The average absolute deviations (in peccent) are: 
DEVMON = 12.9, DEVHYP = 19.9, and DEVGDE = Zara: 

In general, the accuracy of the approximations 
decreased. However, when compared to Table 10, the monova- 


oe 


riable formula improved its performance in 6 out of 15 
cases, and the CDF formula in 7 out of 15 cases, whereas the 


hypergeometric formula lost accuracy in 14 out of 15 cases. 


i“The additional letters in Coltimms DE VMON on ecard 
DEVCDF indicate: 


I - improved 
E - same level 
W - became worse 


4 Q 


it can be concluded that an increase of the jnaumber repair 
Geyeves has a larger improving influense on the accuracy of 
the monovariable andthe CDF formulas than it has on the 
accuracy of the hypergeometric formula. The accuracy of the 
hypergeometric formula iS more sensitive to the ratio of N1 
to N2 than the two others; if performs best when Ni = N2. 

The impact of increaSing n can quickly pe seen when the 
average absolute deviations of Table 11 are compared to 
fame: 9. 


DEVMON DSV AL? DEVCDE 
Pale 9; 7427 60.1 Hye 
Table 11: 239 14.9 Zia 


iaolbe oF dua lleedisplay results for very Similar paran- 
eter values. In corresponding rows all parameters vaiues are 
Hoent ical except Wihyn2 andy. However, in each row of 
Table 11 there are increased but unegual numbers of repair 
cycles, the Y's of Table 11 are larger than in Tabie 9 ia 
Gadel TONeKeCeDmene Latio ¥/MU at about the same levei for 


corresponding rows. 


G. SIXTH SIMULATION SETTING 


Up to this point the timing of demands has been assuied 
Rene at the start or the repair cycle. However, this nay 
not always be the case at the NARF; the timing is dependent 
Gn the length of time for teardown and reassembly of a 
repairable. It is quite possikle that demand for an itena 
could occur at the end of a cycle when the reapairable is in 
its final stage of reassembly. iieadaneror, iLiechapter i I 
we illustrated the impact of the timing of the demands 
during the repair cycle on the shortage costs with a simpie 


example. 
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Another assumption has been that, in tne event of two 
customers Simultaneously demanding tae last unit in inven- 
tory , the supply department would toss a coin to see who 
gets the item. However, it is not unreasonable that the 
supply department might favor the customer with the higher 
expected demand per period. 

In this simulation setting demands are assumed to occur 
at the end of the repair cycle and, in the case of a tie for 
the last unit, the probability that a customer will get it 
is the ratio of his mean demand for the period to the total 
mean demand from all customers. The Simulations were reruns 
of Table 9 with "Stage 1" involving the timing change and 
"Stage 2" involving both the timing change and the new 


PFIOLLLY tude: 
1. Analysis_of the Simulation Results: 


The average absolute deviations (in percent) are: 


DEV MON DEV HYP DEVCDF 
in fable so: 74.7 OUm lige 
For “Stage: Vl: 128 tog 21.5 
For Stage 2: 10.1 gaat , ae 


These results establish with very high significance the 
decisive influence of the timing of the demands on the accu- 
racies of ail three appproximation formulas. They also 
suggest that, 1£ the supply department would consider a 
priority scheme based on expected demands rather than 
tossing a coin for the last unit, the shortage cost approxi- 
mations could be used with more confidence. 

Interestingly, the occurrence of demand within the 
repair cycle does not matter, when the number of repair 
cycles is identical because the cycles are then of identical 
length too. In the case of an egual number of repair 


cycles, the order in which demands occur will not change as 
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Stage i, Seaeie a2, 
Accuracies shifted additionally 
achieved occurrence Changed 
in Table 9 of demand DELOnLeEy rule 

Pao oe dwt ep -CpE ] D=NON D-HYP D=CDr | D=-MON D-HYP D-Cor 
203 230 506 -31 - 40 10 Z =13 O17 
0 = ie -1 16 2 tS 5 8 4 
14 27 Ze en 12 =: ao 2 = 
256 136 oy Ze ire Da 15 alee 64 
BZ 10 a6 4 a id 1 al) 5 
30 14 aig 10 -4 34 5 =o 28 
a ee ae 3 20 26 5 =6 a 
84 62 111 = = 13 14 -4 =15 11 
98 118 151 =a 5) 7) a 45 -4Q od 
245 216 355 =) ~42 -17 a =o -7 
BZ 20 38 = =A 2 =a - 16 -4 
-24 3 a) =a 34 = 5) y 45 510 
17 4 33 =O - 16 8 =G Sates 7 
al —4 -] 11 8 10 7 4 7 
Ze 34 1S) = a =) 3 13 = 3 


Paem Funes 20,000 replications 


long as demands occur at the same point in time within the 


repair cycles for both demand sources. This fact explains 
the relatively good performance of the approximation 
formulas in simulation settings 2 and 4. In the case of an 


unequal humber of repair cycles the order will generally 
change. 
H. SEVENTH SIMULATION SETTING 


As the number of repair cycies increases, the influence 
of the timing of demands within a repair cycle on the accu- 


racy of the approximations will decrease. 
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Simulation settings 1-6 provijed strong evidence that 
the CDF approximation formula iS Clear bye treme once. it 
does not provide sufficient accuracy ~and will | °eee-7ee 7 
from further simulations and analysis. 

The following restrictions for the parameter values 
should ensure an acceptable accuracy for tn2 monovariable 
and the hypergeometric approximation ‘formulas when they are 
applied to the binomial NARF model. The restrictions are 
derived from simulation settings 1-6: 

1. The number of repair cycles for each demand source is 
greater than or equal to 10. 

2. The mean of the sum of both demands per period pius 
1.50 times its standard deviation is an upper bound 
for Y. 398 

3. The ratio of the laryer unit shortage cost to the 


Smaller iS at most three. 


Table 13. represents results of Simulations which were 


run with parameter values within the above restrictions. 


1. Analysis of the Simulation Results: 


The average absolute deviations (in percent) are: 
DEVHGN ="625 and DEVEYPe==5- 7c. The maximal absolute devia- 
tion is 14 percent for the monovariable formula and 


16 percent for the hypergeometric formula. 

Because the fifteen Simulation runs are independent 
observations for DEVMON and DEVHYP, nonparametric statis- 
tical analysis can be performed. Conover ‘Ref. 8 p.118} 
provides a nonparametric method for computing tolerance 
LIL t Ss His method applied to Table 13 establishes that, 


with a confidence level of 95 percent, at least 81 percent 


,'S5iIf Y is _equal to the mean of the sum of demands per 
period plus 1.5 of its standard d2viations, shortages will 
occur 1n about one out_ of 15 periods. This quantification 
1S a rough estimate derived From the standardized normal 
distribution, 1.¢€. P(Z>15sj=02 67 


uy 


2. able 13: _ Runs with Restricted Varying 


Ni N2 PI P2 BUS teeter Le Seel SHeesc DEVHCN DAVHYP 
elo wOe2 068 14.0 220 16 100 200 909 28 il =) 
epost 0.0 1556 2.7 19 300 100 109 19 5 a 
oie 0.7 O56 1757 2.5 2055100 300 3891 53 14 16 
fom ees YO. 5 1350 2.) 17 300 1007510 11 6 11 
eeze0e i O.9 2059 1.7 23 100 300 175 24 3 2 
meee es 022 961 2.6 12 300 100 175 4 3 0 =| 
ieee so 0-1 1is2 1.6 13 100 250 311 20 =.2 Q 
mimes) 0.9 S450 22537556 100 300 113 DD = =e 
melee 2 0.2 4550, 250, 8 100 300 975 16 ae eZ 
Peon. 20-6 1.072265 16 300 100 305 ho = 12 =) 
eeu s ss 0.59 19238 22.4 23 100 300 129 18 = 1 soil 
es Uevs 7 S229 364 37 | 100 300 610 36 = 2 a 
(eens as S299 27 300° 100 711 6 5 7 
Torelc Oe2z 0s Fez 25.3 86 300 100 117 118 =e 0 
figment Os 30065 2560 569 25 100 300 637 ie -4 = 


Peer run: 20,000 replications 


of tne absolute deviations of the monovariabl2 and hypezrgeo- 
metric formulas will not exceed the maximal absolute devia- 


Evens of Table 13. 


I. SUMMARY AND CONCLUSIONS 


This chapter has shown that the monovariable and tke 
hypergeometric formulas provide significantly more accurate 


Besutes than the CDF ftormula does. The following additional 


16STD denotes the standard deviation for the sum of Dbotn 
demands per period. Sif Te cwOteommtme oa vorage Valuc TOF 
the shortage costs per period as observed during simulation. 
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remarks refer to the monovariable andthe hypergeometric 
£Ormu lass 

All simulations show a ciearly recognizable trend that 
the accuracy decreases as Y increases with levels of Y whach 
correspond to most "real world” situations (12e2))meaa to 
mean +30%). This suggests that the approximations are not 
good in the upper tails of the demand distributions. 

The simulations suggested that the approximations 
provide their most accurate results when there is an iden- 
tical number of repair cycles. Even in the case of a high 
vaCiability of the unit shortage costs (i.e. 3 1,000 and 
$100 ), tke average absolute deviation from the simulation 
results ( simulation settings 2 and 4) is about 10 & for the 
monovariable formula and 5 % for the hypergeometric formula. 

The variability of the unit shortage costs is an impor- 
tant factor for the accuracy of the approximation formulas . 
Reducing the variability of the unit shortage costs’ fron 
$1,000 and $ 100 to $ 200 and 3 100 almost doubies tne 
accuracy to be expected from both approximations, as was 
demonstrated in simulation setting 1. 

Simulation settings 4 and 6 suggest that the approxima- 
tions perforn Significantly better when the demands are 
assumed to occur at the end of the repair cycle. The influ- 
ence of the timing of the demands on the accuracy of the 
approximation formulas decreases as the number of repair 
cycles increases. 

In settings with different numbers of repair cycles per 
period the monovariable and the hypergeometric formulas 
provide about the same level of accuracy. By imposing the 
following restrictions for most practical applications an 
acceptable level of accuracy can be expected: 

1. The number of repair cycles for each demand source is 


greater than of equal ton10- 


¢ 
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2. The mean of the sum of both demands per period plus 
1.50 times its standard deviation is an upper bound 
LOE €ie 

3. The ratio of the larger unit shortage cost to the 

smaller is at nost three. 
More specifically, when observing these restrictions it 
is to be expected that, With a confidence level of 95 
percent, at least 81 percent of the actual absolute devia- 
tions of the monovariable formula will not exceed 14 
percent, and, in the case of the hypergeometric formuia, 


they will not exceed 16 percent. 
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VI. MORE THAN TWO BINOMIAL DEMAND SOURCES 


A. PRELIMINARY REMARKS 


This chapter analyzes the situation when there are more 
than two demand sources and presents results from the asso- 
ciated simulation runs. 1D pantucular, it Seonsiasus sence 
approximation of the distribution of ~—fne SMES e seen anaoee- = 
period by the normal distribution. 

The computational effort required to use the hypergeone- 
tric formula increases very rapidly with an increasing 
number of demand sources. Therefore only the nonovariable 
formula will be evaluated for the case of more than two 


individual demand sources. 


B. METHODS FOR ESTIMATING THE EXPECTED SHORTAGE COSTS PER 
PERIOD 


Two basic ways for estimating the expected snortage 
costs per period in the case of "many" binomial distributed 
demand sources are analyzed in this section. The first is 
to simulate the complete period for a sufficient number of 
replications to provide the desired precision. fe the 
desired precision has been prespecified, the number o£ 
necessary replications can be determined by the progran 
during its execution. This simulating approacn gives an 
unbiased estimate for the expected shortage costs per period 
and, additionally, allows for varyiny the time at which the 
demands are assumed to occur within the repair cycle. 

The second way to estimate the expected shortage costs 
pec period is to apply the monovariable approximation 
LOU MUclar However, before applying the monovariable formula 
the distribution for the aggregate demand has to be 


computed. 


48 


in 


Ene Case of ‘“many™ “binotially distributed  deaanc 


sources with different p parameters there are essentially 


four approaches for computing the distribution of aggregate 


demand: 
1. 


the sum of demands per period by simulation. The 
demand occurrences during the period are simulated 
and added up to yet one observation for the sum of 
demands per period. This is repeated for a specified 
mamber of replications. The observed empirical 
dastribution for the sum of demands during simulation 
can be taken as an eStimate for the distribution of 
the sum of demands per period. 

The second approach is to compute the ex 
tion of the sum of demands. However, 
computations for computing the exact distribution 
grows as a factorial expression whereas the number of 
computations performed during simulation only grows 
approximately linearly along with “the number of 
demand sources. For instance, if there are five 
demand distributions with ¢6,7,8,9,10> repdair cycles 
per period, then a five-fold convolution is reyguired 
which incorporates 7x8x9x10x11 = 55,440 iterations. 
There are two analytical approximating nethods which 
can reduce the huge number of computations invoived 
in finding the distribution of the sum of the demands 
per period. One approach is to approximate the indi- 
vidual binomial distributions by Poisson 


Ge Stributions. 


The second analytical method is to approximate the 
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C. APPROXIMATION OF THE SUM OF DEMANDS BY THE NORMAL 
DISTRIBUTION 


An advantage of the monovariable approximation formula 
is that it does not reguire the consideration of ali 
possible permutations. It requires only the iistribution of 
the sum of demands per period. Therefore it can be applied 
when the individual binomial distributions are approximated 
by normal distributions. The sum of normal random variables 
is again normally distributed. Its mean and variance are 
the sum of the individual means and variances [Ref. 9 
p.- 405 }. 

The Poisson distribution could also be used as an 
approximation since the sum of Poisson distributed randon 
variables is a Poisson random variable whose wntean 1s the sun 
of the individual means. [Ref. 9 p.405]. 

The advantage of the normal distribution over the 
Poisson is that it has two parameters. One parameter deter- 
mines the level (=mean) and the otker one the _ spread 
(=variance). The effects of the two parameters are putually 
independent and permit the modelling of a wider variety of 
demand behavior. Therefore the normal distribution is often 
preferred over the Poisson distribution when one attempts to 
fit a theoretical distribution to an empirical distribution. 

A continuity correction was used to reduce the error of 
approximating a discrete distribution with a normal distri- 
bution. If, for instance, the probability associated with 
the integer I is to be computed, the normal approximation 
with continuity correction uses the area under the normal 
density from (I - 0.5) tO - (tet 0 5) to approximate the 
proba basi. 

A further reason for using the normal distribution is a 
conseguence of the Central Limit Theoren. Prom Ene Gen tear 


Limit Theorem it follows that as the number of demand 
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sources increases, the aggregate distribution of demands per 
period will be increasingly better approximated by the 


POEMal distribution. 


De. RUNNING THE SIMULATION 


In case of an arbitrary number of demand distributions 
the monovariable approximation formula useS an approximating 
distribution to evaluate the expected number of shortages 
per period. Therefore, deviations between the result of 
Simulation and the monovariabie approximation formula when 
aggregate demand 1S assumed to be normally distributed might 
be caused by both the process of applying the approximation 
formula and this use of the normal distribution. 

The output tableau of the simulation program is hbasi- 
cally an extension of the output tableau described in 


Chapter 4, section D of this thesis. To enable an efficient 


Output analysis, the program provides more key statistics 
eh nD the Former PEOgLan Lor two binomial demand 
distributions. 

NOFDMD= 21 #OF TIME UNITS= 302 REPLI=50000 I-SEED= S09 
AOE REP PROB A- UN DiS Shon Ts Eile 
GYCLES BELETY MEAN EST N¥P N*¥Q 

1: 10 lem 0 20 Ue 2200 eo. OW De Oy cs 0 

ue (lS 2 eo 22 ec 0 0 ee I 0 ear 30 

TOTAL MEAN= 14.00 TOIAL gs DEV= 2.008e) Y= 156 MAX-DMD= 25 
SIMULATION: AVE-COST= Zawso- Or SHOG'S= 0215 25.¢ 

SHoriu= W5307 SD-COSi- Oe eed FOr —sH'S= > 0.00237 

NORMAL-APPROXIM.: AVE-COST= SOON #OF SHOR'S= 0.1616 

SUMMA to rALISTICS: 

SD/AVE= 1.47%8 PISIMA= 185.9759 PINORA= 135.715 %t° 

DEVMON= 5-83%11 DEVNSH= 5.97412 DEVPIA= -0.14%15 
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Number of demand Sources; 

Number of time units per period; 

The last two columns indicate how well the individuai 
binomial distribution might be approximated by an 
normal distembucion In jeneral, nape andmen= (1—p) 
should both have at least a value of five to ensure a 
sufficient fit Sion the individual distribution. 
However, applying the monovariable formula incorpo- 
rates only the distribution of aggregate demands. 
With an increasing number of demand sources it is to 
be expected that according to the Central Limit 
Theorem the approximation of the aggregate denand 
distribution will be by far better than indicated by 
the N*P and N*(1-P) values for the individual 
distributions. 

Standard deviation of the sum of demands; 

Estimated expected shortage costs per period; 

Average number of shortages per period as they were 
observed during Simulation. 

The observed standard deviations for the estimates 
give a numerical hint about the precision of the 
estimates. 

The observed standard deviation of the sample mean 
divided by the sample mean of the shortage costs per 
period. 

Average unit shortage costs derived fron the simula- 
eiLOl i.e. accumulated observed shortage costs 
divided by the observed number of shortages. 

Average unit shortage costs weighted by the means of 
the individual demand distributions corresponding to 
the way the monovariable approximation formula uses 
these weights. 

Deviation of the result of the monovariabie approxi- 
Mation formula from the estimate for the expected 


shortage costs per period derived from simulation. 
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12. The deviation of the average number of shortages used 
by the monovariable approximation formula from the 
average humber of shortages aS observed during 
Simulation. 

13. The deviation of the mean weighted average unit 
shortage cost computed when applying th2 monovariable 
approximation formula from the observed average unit 
shortage cost. VE vist and DEVPIA provide a rough 
means for distinguishing whether the deviation 
between simulation results and monovariabie approxi- 
Mation is caused by the approximation formula or by 
the approximation of the distribution of the aggre- 


gate demands per period by the normal distribution. 


E. SIMULATION RESULTS 


This section presents conclusions derived fron sinmula- 
tion runs performed in order to evaluate the performance of 
the mohnovariable approximation formula when the total denand 
distribution is approximated by the normai. For all simula- 
tions presented in this section the range of possibie paran- 
eter values is within the limits of the three restrictions 
listed at the end of Chapter V. 

Table 14 is a direct continuation of Table 13. The only 
difference is that the distribution of the sum of demands is 
how approximated by a normal distribution. 

The average absolute deviations (in percent) are: 
DEVMON = 11.1, DEVNSH = 10.4, and DEVPIA = 5.8. The ‘DEVPIA' 
column of Table 14 is almost iijentical to the ‘'DEVMON' 
column of Table 13. The ‘DEVNSH' column of Table 14 measures 
the additional inaccuracy incurred by applying the noraal 
appEOxINdtwon instead of using the exact distribution for 


the sum of demands as was done in Table 13. 
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1. Table 14: _ Runs with Identical Parameters as 
Table 13 '’ 

Ni N2 Pl mee MO STD Y Pail Pi2 Seed DE VMole DEV NSiheore ye a 
10 15 0.2 0.8 14.0 220° lomo ee 2 Coe Diese 6-50. Oem 
12°18 0.4 026 15°26 227 IoC eee Ome 9.4 5.4 328 
15 12 0.7 O26 17.7 2.25520) | 0Gmciac?) 22a WOES A a6 
10 20 0.3.0.5 13.08257 17) CU ieee eee Sal alien Gag 
11 22 081 0.9 20.S9mie7 23 See Oee ete Ps ONS Del eo olan 
271 14 0.4.3. 052 921 226612 SCC © ao ioull St 
12 16 0.8 0.1 Ada2 VeCw ee Cum oe call Zo  1.Oeno 
20 40 0.9 0.9 56.0 22.356) 1COme Oe =a Nile 2 
10 15 0.2 0.2 5,0 2.0 8 100° 300° 975 9-23.) eee 
20 15 0.2 0.6 13.0 2.6 16 300 1007305 5s = eae 
12.18 0.9°055 19506254 23 OCR CCS = ae 19 6:4 aoe 
13 39 0.4 0.7 32.5 3245527 eee Oe Se. Cre Geter 
11 33 .05 0.1 3.9 1.9 7 300 1000011 —333 cee 4.2 
19.12 0.2 0238 7.2 2o55 Roe OC tar ee eo om 
10 40 0.3 0.5 23.0 3.35925) 100) 0c = Seng = 1. ee cee 
EaGh run: 20,000 replicarvoms 

If in one row of Table 14% the actual values of DEVNSH 
and DEVPIA are either both negative or both positive then 


the actual overall deviation (DEVHON) 
DEVNSH and DEVPIA. Otherwise, when tne Signs of the 
DEVNSH and DEVPIA are different tnen the 
aggregate affect of DEVNSH and DEVPIA on the value of DEVMON 


results 


represents an accunula- 
Lomo 
actual values for 
in a partial compensation of their individual 
In Table 14 the actual EOE SOE Vion 


DEVPIA have different signs in nire out of fifteen rows. 


affects. values and 


17DEVNSH 


= deviation in number of Short'"s(in percent) ; 
DEVPIA = 


deviation in average unit shortage cost 
(2n percent). 


54 


The average absolute deviations of DEVNSH and DEVPIA 
indicate that about two thirds of the overall average abso- 
lute deviation is caused by using the normal distribution as 
an approximation for the distribution of the sum of demands 
and the rest from applying the monovariable approximation 
formula. It can be concluded that in the case of two demand 
sources it 1S worthwhile to compute the exact distribution 
of the sum of demands if the monovariable approximation 
formula iS going to be applied. 

Simulation results for three denmard sources are 
displayed in Appendix C. The mean of the distribution of 
the sum of demands is equal to 21.6 demands per period and 
its standard deviation is egual to 3.75. Tnese tables show 
that the accuracy of the overall approximation decreases 
with increasing Y. The following results (in percentages) 


were obtained: 


DEVMON DEVNSH DEVPIA 
= 20 a is. Om = 1) ot) 
= 23 =-3.4 eZ mee 2 
= 26 -7.4 Le: roe! 


These results indicate again that the accuracy decreases 
as Y increases. 

Appendix D contains results of fifteen simulations for 
three demand sources. The fifteen cases of Table 14 were 
expanded by one additionai demand source in orier to provide 
a basis for comparing the results of Appendix D with Table 
14. For each tableau of Appendix D the ratio Y/M4U is about 


the Same aS in corresponding rows of Table 14. 


ie average absolute deviations (in percent) are: 
DEVMON = 5.8, DEVNSH = 4.2, DEVPIA = 3.8. For Table 14 they 
were: 11.1, 10.4, and 5.8, respectively. Th2 expansion to 


three demands sources shows already a strong numerical 
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influence from the Central Limit Theorem. The average abso- 
lute deviation for the number of shortages dropped from 
1024. tomes 27... In six out of fifteen cases tne actual 
values of DEVNSH and DEVPIA have different signs which 
reduces the overall deviation for those cases. 

Conover'ts method [Ref. 8 p.118] for computing tolerance 
limits applied to Appendix D establishes that, with a confi- 
dence level of 95 percent, at least 81 percent of the abso- 
lute deviations(DEVMON) will not exceed 11.9 percent if the 
three restrictions for the ranges of parameter values are 
applied. This accuracy should be acceptable for most 
applications. 

Appendix E contains results of simulation runs for five 
demand sources with one parameter set and three levels of Y, 
i.e. mean, mean plus one standard deviation, and mean plus 


1.5 standard deviations. The deviations (in percent) are: 


DEVMON DEVNSH DEVPIA 
= 42 oe Ue0 va 
= 46 0.5 ie ley 
= 48 Ores Se PO Lae 


These results indicate that an acceptable level of accu- 
racy can be expected in the cases of five and more demand 


sources. 


Fe. SUMMARY AND CONCLUSIONS 


This chapter has shown that in case of two demand 
sources it 1S worthwhile to compute the exact distribution 
of the aggregate demands. In the case of three demand 
sources, the benefit gained from computing the exact distri- 
bution for the sum of demands is significantly lower. 
However, if the numbers of repair cycles are small - which 


means a relatively small computational effort is reguired in 
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computing the exact distribution of the sum of demands - it 
is recommended that the exact distrinution for the sum o£ 
demands be computed instead of using the tnormai approxima- 
tL Ohl Otherwise for the cases of three or more demand 
sources the normal approximation should be applied. 

All simulation results presented in this chapter indi- 
cate that the sign of the actual deviations caused fron 
approximating the distribution of the aggregate demands 
(DEVNSH) 1s independent from the sign of the deviation 
resulting from applying the monovariable approximation 
formula (DEVPIA). If for a given parameter setting the 
Signs of DEVNSH and DEVPIA are different then the overall 
deviation (DEVMON) represents a partial compensation of the 
aggregated affects of DEVNSH and DEVPIA; otherwise, the 
overall deviation is about egual to the sum of DEVNSH and 
Dave LA. Therefore, the impact of the normal approximation 
on the variability of the accuracy will be stronger than its 


impact on the expected level of accuracy. 
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VII. SUMMARY AND SUGGESTIONS FOR FURTHER ANALYSIS 


A. SUMMARY 


Chapter II proved that the monovariable approximation 
formula for the expected shortage costs per period provides 
the exact value for the expected shortage costs per period 
in the case of any number of Poisson distributed individual 
demands. Because the monovariable approximation does not 
provide the exact value for the expected shortage costs per 
period in the case of binomially distributed individual 
demands as they are modelled by tn2 NARF model, two further 
approximation formulas, the hypergeometric formula and the 
CDF formula, were introduced in Chapter III. 

Chapter IV described the simulation methodology which 
was used in Chapter V to evaluate the performances of the 
three approximation formulas in the case of two binonmially 
distributed demand sources. Chapter V showed that the mono- 
vaciable and the hypergeometric formulas provide about the 
same level of accuracy. Both formulas provide significantly 
more accurate results than the CDF formula. They are there- 
fore the ones that should be used in applications. 

Chapter V showed that three factors are significant in 
cetermining the accuracy of the three approximation 
fOr muvac- 

1. The occurrence of demand at the end ofr the repair 
cycle improves the accuracy of the approximations. 
However, in almost all "real world" situations which 
are modelled by the NARF model, demand will occur at 
the bejinning of the repair cycle. The approximation 
error incurred hereby will decrease as the numbers of 


repair cycles increase. 
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2. The accuracy of the approximations for the expected 
Shortage costs per period decreases as the available 
Guantity Y increases reflecting the fact that the 
approximations are not so good in the upper tails of 
Bier arstributaons. 

3. If both unit shortage costs were the same, the 
problem of approximating the expected shortage costs 
per period is reduced to approximating the expected 
number of shortages per period. The higher the vari- 
ability of the unit shortage costs, the less accurate 
1s the approximation for the expected shortage costs 
per period. 

By imposing three restrictions on the range of parameter 
values, an acceptable expected level of accuracy can be 
achieved for the monovariable and the hypergeometric formula 
in the cases of two demand sources. These restrictions are: 

1. The number of repair cycles for each demand source is 

greater than or equal to 10. 
ze The mean of the sum of both demands per period plus 

1.50 times its standard deviation should be an upper 

bound for Y. 18 

3. The ratio of the larger unit shortage cost to the 
smaller is at most three. 

Tolerance limits for the accuracy were computed frou 
Simulation results in Chapter Vs. With a confidence level of 
95 percent, at least 81 percent of the actual results fron 
the monovariable approximation formula will incur at most an 
absolute deviation of 14 percent from the exact numerical 
value of the expected shortage costs per period, and at most 


16 percent in case of the hypergeometric formula. As a 


ee  —— —————ee 


2 tec gialwtoO the mean sof the sumeot demands per 
period plus 1.5 of its standard deviations shortages wili 
occur in about one out. of 15 periods. ThES juagwti fication 
is a rough estimate derived from the standaftdized normal 
SrotrabwrlOn, 1.2ec. P (Z>1.5)=0. 067. 
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general tendency, the hypergeometric formula “seens ico 
perform slightly better than the monovariable approximation 
formula, especially in the case when there are equal numbers 
of repair cycles associated with each demand source. 

Chapter VI expanded the setting to more than two binomi- 
ally distributed demand sources. The problem with the 
hypergeometric formula is that all possible permutations of 
demands must be evaluated whereas the monovariable formula 
reguires only the distribution of the agyreyate demands per 
period. Because the number of computations needed to apply 
the hypergeometric formula increases rapidly as the number 
of demand sources is increased, only the monovariable 
formula was evaluated in the cases of more than two demand 
sources. 

Chapter VI showed that in case of two demand sources it 
is worthwhile to compute the exact distribution of the 
agjregate demands rather than using the normal distribution 
to approximate the distribution of agjregate demands. in 
the case of three demand sources, the benefit gained fron 
computing the exact distribution for the sum of demands was 
Significantly lower .due to the increased computationai 
effort required. However, i1f the numbers of cepair cycles 
are small, that is not much of a problem and it is recon- 
mended that in the cases of three or four demand sources 
that the exact distribution for the sum of demands be 
computed. Otherwise for the cases of three or more demand 
sources the normal approximation can successfully be applied 
if certain parameters are held within a limited range of 
values. 

All simulation results presented in Chapter VI indicated 
that the signs of the actual deviations caused by approxi- 
Mating the distribution of the aggregate demands were inde- 
pendent from the signs of the deviations resulting fron 


applying the monovariable approximation formula. Try eons 
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one specific parameter setting, the signs are different then 
their effects on the overall deviation get partially compen- 
sated. Otherwise, the overall deviation is about equal to 


the sum of the two single deviations. 


B. FACTORS TO BE CONSIDERED FOR PRACTICAL APPLICATIONS 


In case 9f binomially distributed demand sources’ the 
following factors should be considered. 

The NARF model incorporates mutually indepandent binoni- 
ally distributed demands. AS waS indicated in Chapter II, 
the Poisson distribution can be used to approximate the 
binomial distribution. The approximation provides an accep- 
table level Of aceumdacy ini and | ps0. 05; (Cee BS 
N2100>5and nps10 [Ref. 10 p.79}. If the actual demand 
distributions at the NARF all exhibit such parameter values 
then the monovariable formula can be used along with the 
Poisson distribution for the aggregate demanis. 

AS alast resort, Simulation of the actual demand 
behavior may be the oniy approach for oaetermining the 
expected shortage costs per period. Mimearticular, 1: one 
of the restrictions listed above is violated, simulation 


should be used to estimate the expected shortagjye costs. 


C. SUGGESTIONS FOR FURTHER ANALYSIS 


The approximation for the expected shortage costs 
appeared to be more accurate when the timing of demand noved 
to the end of the repair cycle and when the priority disci- 
pline was based on expected demand. in the case of the 
latter result, a minor change in the policy of running the 
system is suggested. 

As the stock gets more and more depleted during the 
period a rational policy would be to reserve stock for 


customers with the higher unit shortage costs and expected 
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higher demands. The higher the variability of the custoner 
specific unit shortage costs and the lower the holding cost 
pec stocked item, the greater is the benefit from switching 
from a "first-come, first-servea" rule to a fore sophisti- 
cated ruie which allows for the different economic impor- 
tance of the different demand sources. Analysis has to show 
whether the costs for implementing this flexible policy are 
justified by cost reductions in the sum of economic costs 
associated with running the inventory systen. 

In the case of the former result, further study of the 
influence of timing 1S warranted. Perhaps a shortage cost 
approximation can be developed which includes the timing as 


a parameter. 
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APPENDIX A 


== ce ce ee 


COMPUTER PROGRAM FOR SIMULATING TWO BINOMIAL DEMAND SOURCES 


AANA A AAA AA AAAA ANA ANAANAA AAA AAA AAA AAANAAAANANAAAANNAANANAAAAAANNAIANANM 


NAME: WILHELM F. NOEGGERATH 
DATE aewUbayeionn 1984 
BINOMIAL DISTRIBUTION, TWO DEMAND SOURCES 
- SIMULATION OF BERNOULLI TRIALS PER REPAIR CYCLE 
- MONOVARIABLE APPROXIMATING FORMULA 
Tele eRrGEeCieteec OhvenrArtsD APPROXIMATION FORMULA 
=—"CUF OQORLENIAZED APPROXIMATION FORMULA 
VARIABLES: 
APPROX RESULT OF MONOVARIABLE APPROXIMATING FORMULA 
AUX... AUXILIARY VARIABLES 
AVECOS AVERAGE COST PSR PERIOD 
Die Lie ECG Ea teoecLAS LF TWO DEMANDS EN ONE TIME UNIT 
AND UGHNE, UNIT. OF STOCK LEFT ONLY 
COsCpr RESUPT VOL GDr—- ORTENTATED APPROXIMATION FORMULA 
SUsSeri, COST PURmeeekENT STIMULATED PERIOD 
COseuUn CUD ULA Tepe eCOST FOR ALL SIMULATED PERIODS 
COSutP RESUC Teer enyYPERGEOM. ORIENTATED APPROXTMS FORMULA 
Coss” ) Culver PERIOD SOUARED COSTS 
ek COUNTS PERIODS IN WHICH SHORTAGES OCCUR DURING 
SIMULATION 
Diy CPE SD Ev Arle Ne IN PZRCENT OF RESULT CDF ORIENTATED 
APPROXIMATING FORMULA FROM RESULT OF SIMULATION 
DEVEYP DEV. IN PERCENT RESULT HYPSEO. APPROXIM. FORMULA 
DEVHOH DEV. IN PERCENT RESULT MONOVAR. APPROXIM. FORMULA 
DMD SCALAR, CUMULATIVE DEMAND IN ONE PERIOD 
DMD1 ARRAY OF DEMANDS WITHIN ONE PERIOD CAUSED 
BY END ITEM #1 
DMD2 ARRAY OF DEMANDS WITHIN ONG PERIOD CAUSED 
BY on LT Bia Z 
itv eiGniMOn eng enin CYCLE SPOR TEND ITTEMN i] 
miylviG2 LENG? HVOr REPAIR CYCLE FOR END ITEM 42 
Psp ONL EALS SEED, GETS PRENTED 
PICT ely erties TOR AHiCH END ITEM INS CURRENT SIMULATED 
PERDOD IASI DEMAND SATISFIED 
MU Sime, Ul + MUZ 
MU1 MEAN BINOMIAL-DIST. DEMAND CAUSED BY END ITEM #1 
MU2 MEAN BINOMIAL-DIST. DEMAND CAUSED BY END ITEM #2 
N NUMeR or DISCRETE PIIMs UNITS PER PERTOD 
Pio awes N MUST PE A MULTIPLES OF NI AND OF NZ eee st 
N1 NUABER OF REPAIR CYCLES END ITEM #1 IN ONE PERIOD 
N2 NUMBER CF REP ADRe! CLES SENOe ITEM #2 IN ONE PERIOD 
Pil SEORTAGE COST ITEM #1 
PE SHORTAGE  C@sl LLEM YZ 
PIE MNO Deno o FCl FOR DUD) PER PRRIODSEGAUS ED 
Di cel iets # | oAND #2 
PMF1 PROBABILITY MASS PCI SEGR DND PER PERIOD 
CAC or Dest Ee Maral 
Pir Z PROBABLUIAY NASperel FOR DID PER PERL OD 
Cause D BY. E=-LTEM FZ 
P1 PROBAB. OF DEMAND CAUSED BY END ITEM #1 
PN ONE ene Pah CY el r 
P2 PresAbewOr DEMAND CAUSED BY END ieee +t 2 
PONE whe PALR CYCLE 
Rel DieenuaBER OF SIMULATED PERTODS 
SALOND SATISFIED DEMAND 
SDCOS ESTIMATED STANDARD DEVIATION OF AVERAGE COST 
Peer oh ea) 
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) AND. (Y.GT.0)) 
: 


Y 
GSCUMTCOSEHT 


SATDMD 
J + 
J 
C 
GO ssOet CO Osen iE. «2 


ENDWHILE 


ENDIF 


J 

COSCUM = 

COSS® = 
CONTINUE 
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HU1 = BISELOAT (NI) 
MU2 = P2*FLOAT {N2 
MU = MU1 + MU2 
COMPUTATION OF SAMPLE MEAN AS 4EAN ESTIMATE AND 
COMPUTATION OF ESTIMATED STANDARD DEVIATION OF 
MEAN ESTIMATE 
AVECOS = COSCUM/FLOAT (REPLI 
AUX = COSCUM*¥*2/FLOAT (REPLI 
AUX1 = FLOAT REPLI) #FLOAT (REPLI-1) 
SDCOS=S9RT ((COSSQ-AUX) /AU 
,FRITE (9,601) REPLI, N, ISEED, N1,P1,MU1,PI1, : P2, 
gf 
601 FORMAT (" 71(**') 7,3 REPLI ATIONS: 
* TIME UNITS PER ERIOD:? tint it 8 fofe 7 TE 
os OF KEPCYCLS_1:',13 ees MU1=t ea ; 
*  F7.2,T67,! =" 147/, UN@On BAPGNGrSe7 wey 
* ot  po=t Pelu, MU 2= eee eee PI2='(F7.2,/ ) 
WRITE (9-602) AVECOS, SDCOS 
602 FORMAT('" ',"RESULT OF BERNOULLI TRIALS FOR EACH!, 
x ' REPAIR CYCLE:'/, 
+ ' ESTIMATED MEAN OF SHORTAGE COSTS PER', 
x ' PERIOD', F8.3,/ 
*' ESTIMATED ST-DEV.’ FOR ESTIMATE OF THE MEAN:',F8. 3) 
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APPLYING THE MONOVARIABLE APPRKOXIMATING FORMULA AND 
EVeeUnTONGeLHE EXACT EXPECTED FRACTION OF PERIODS 
ee cs SHORTAGES OCCUR 
Whe E { die Les NI i 
PMFT(I+1) =DGA uA (DPLOAT (N14 1)) / 
* NEMEC Tey +1)) * 
* See fay i) 
MF1(I+1) = PMF1 T+) *(DBLE (b1)) **1 * 
(DBLE(1.0 - P1))**(N1-1) 
ne I+1 
40s 
WHILE a Ee N 2 
PMF (I+ ) =DGA MRC EC OUn Oe 
Wepgutua or oa (N2-I+1)) * 

+ GAMMA DELOAT (1+ 1))) 
PMF2(I+1) = PMF2(I+1) *(DBLE(P Daa = 
ie— +4 

ENDWHILE 

I = 0 

Wee ieeeeeoN Lt 2.) 
PMF(I+1) = 0.0D0 
I=1I+ 7 

ENDWHILE 

I = 0 

WHILE (2 Seen} 
WHILE ¢ LBs L 

BMF 19+ 1) ee +J+1)+ PMF1(I+1) * PMF2(J+1) 

ENDWHILE 
I=1I+ 1 

ENDWHILE 

SUM = 9.0 

SHEXAC = 0.0D0 
Y s+ 1 

WHILE (I. LE.N1 + N2 
SUM = Sut + FLOAT (I- Y)#SNSL (PAF (I + 1)) 
SHEXAC = SHEXAC + PMF (I+1) 

Peet 4 

ENDWHILS 


APPROX= eee (02) Wea teH02}) +508 


604) APP 
7 RESULT OF MONOVAR. * 
*hDPROXIMATING FORMULA:',F8.3,/) 
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606 


CALL HYPGEO (N1,N2, PHF 1, PMF2, Pll, Pi2 eee 
CALL SUBCDF (N1,N2, PHF I) PMF2) Pil Pcie ee) 


COMPUTATION OF SUMMARY Stes 


SHEXAC = SHEXAC*100.0 
SHSIMU = (CTR/PLOAT (REPLI) ) *100.0 
IF (AVECOS.GT.0.005) 
DEVCDF= ( COSCDF /AVECOS - 1 - 9) 2190.0 
DEVHYP= ¢ CCSHYP /AVECOS - 1.0) *100.0 
DEVMON= (APPROX/AVECOS - 1.0) *100.0 
SDMEAN = {SDCOS7AV ECOS) #190. 
WRITE (9,605) MU,SHEXAC,SHSIMU,SDMEAN, DEVMON, 
x DEVHYP,DEVCDF 
FORMAT (* SUMMARY: ere - 
x F6.2,°% SHSIMU=',F6. 2 "KY SL 
x ' SD/MEAN=',F6.1,'% DEVAoN= (Pb. 1, 
ee 39,"DEVHYP=',F6.1,1753,'%  DEVCDF=',F tg ) 
WRITE (9,696) MU, SHEXAC, SHSIMU 
FORMAT (* SUMMARY: He Wane ee SHEXAC=', 
x F6.2,'% SHSIMU=" FG. 2," 5 ee 
* ' SD/MEANS 3% Ex DEVMON= ***x4keeet 733, 
* 'DEVHYP= **¥*xexex' TS] 'DEVCDF= sok sok TY 
ENDIF 
STOP 
END 
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RETURN 
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RETURN 
END 


COMPUTER PROGRAM FOR SIMULATING N BINOMIAL DEMAND SOURCES 


$JOB 

c NAME: WILHELM F. NOEGGERATH 
C DATE: JULY, 21ST 1984 
e BINOMIAL DISTRIBUTION, N DEMAND SOURCES 

C - SIMULATICN OF BERNOULLI TRIALS PER REPAIR CYCLE 

e — APPLYING MONOVARIABLE APPROXIMATION FORMULA 

C BY USING NORMAL APPROXIMATION TO BINOMIAL DIST'S 
C VARIABLES: 

Gueee AUX. AUXDLUARY VYARTASLES 

C AVECOS AVERAGE COST PER PERIOD 

C  COSCRT COST FOR CURRENT SIMULATED PERIOD 

C COSCUM CUMULATED COST FOR ALL SIMULATED PERIODS 

C COSSQ CUMULATED PERIOD SQUARED COSTS 

C CRIDMD VECTOR WHICH STORES ORIZINAL DEMAND OF CURRENT 

C TIME UNIT 

C CUMDMD CUMULATES DEMAND DURING PERIOD 

C CTR COUNTS PERIODS IN WHICH SHORTAGES OCCUR DURING 

c SIMULATION 

C  DEVMON DEVIATION OF RESULT OF MONOVARIABLE APPROXIMATION 
C FOR EXPECTED SHORTAGE COST PEFR PERIOD FROM 

€ APPROXIMATION PESULT IN PERCENT 

C  DEVNSH DEVIATION OF RESULT OF MONOVARIABLE APPROXIMATION 
c FOR EXPECTED NUMBER OF SHORTAGES PFR PERIOD FROM 

é SIMULATION RESULT IN PERCENT 

C  DEVPIA DEVIATION OF THE AVERAGE OF THE UNIT SHOETAGE 

é COST AS USED IN THE NORMAL APPROXIMATION FROM 

ec THE AVERAGE UNIT SHORTAGE COST AS OBSERVED 

é DURING SIMULATION 

C DMD DEMAND MATRIX FOR ONE PERIOD 

e ROW=DMD-SOURCE COL=TIME- UNIT) 

C INTVL VECTOR WITH LENGTH OF REPAIR CYCLES FOR 

C DIFFERENT DMD SOURCES 

G IN NUMBER OF DISCRETE TIME UNITS PER PERIOD 

C *k*x N MUST BE A MULTIPLE OF EACH ELEMENT OF REPCYC *** 
C NOFDMD NOMBER OF DIFFERENT DEMAND SOURCES 

C NOFSHN AVERAGE NUMBER OF SHORTAGES PER PERIOD 

C DERIVED FROM NORMAL APPROXIMATION 

C NOFSHS AVERAGE NUMBER OF SHORTAGES PER PERIOD 

e CBSERVED DURING SIMULATION 

a Ff VECTOR , PROBABILITIES OF DEMAND 

C IN ONE REPAIR CYCLE 

Ge pa VECTOR . UNIT SHORTAGE COSTS FOR DIFFERENT 

G DEMAND SOURCES 

C  PINORA MEAN WEIGHTED AVERAGE UNIT SHORTAGE COST 

C  PISIMA AVERAGE UNIT SHORTAGE COST AS OBSERVED DURING 

C SIMULATION 

C  REPCYC VECTOR WITH #0F REPAIK CYCLES PER PERIOD 

C REPLI NUMBER OF SIMULATED PERIODS 

C  SATDMD CUMULATES SATISFIED DEMAND FOR CURRENT PERIOD 

C  SDAVE ESTIMATED STANDARD DEVIATION FOR THE MEAN 

G ESTIMATE FOR THE EXPECTED SHORTAGE COST PER 

é PERIOD DERIVED FROM SIMOLATION DIVIDED BY THE 

é ESTIMATE I.E. DIVIDED BY THE SAMPLE MEAN 
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SIMULATION RESULTS FOR THREE DEMAND SOURCES, 


Y INCREASING 


NOFDMD= 3 #OF TIME UNITS= 72 REPLI=20000 I-SEED= 901941 
#OF REP- PROBA- UNI IT—SiGr. EEL 
CYVCinPees BL Ldely MEAN COS N* Pp N*QO 

1: 18 1: O20 1: 3.60 1: 300.00 3 Some ed 
2: 36 2: 0.40 2: 14.40 2: 300.90 loo? 1.6 
33 le 3° 50230 Be S00 3: 100.00 Beo o.4 

TOTAL MEAN= 21.60 TOTAL S-DEV= 3.75 Y= 20 HAAX-DMD= 66 

IMULATION: AV E-COST= 6 55a SOF. SHO k Vase 2<4178 

SHSIMU= 61.06% SD-COST= Se 27 SD-#OF SH'S= O70 195 

NORMAL-~APPROXIM.: AVE-COST= 646.67 #OF SHOR'S= 224259 

SUMMARYS STATIS TACs: “ 

SD/AVE= 0.80% PISIMA= 2102955 PINORA= 256.678 

DEVMON= -1.29% DEVNSH= 0.303 DEVPIA= -1.58% 

NOFDMD= 3 #OF TIME UNITS= 72 REPLI=20000 I-SEED= 901941. 
#OF REP- PROBA- UN tt— SHOR: Fit 
CY Ciges Ba bey MEAN Coo. N*¥P WW * O 

loom ie ce Ore O es 3200 le 300.00 320, 14.48 
2:6 ee ean Ze 14.40 Ze 300.09) 1Q24e 2 126 
320 te Soe as 3500 SP 100.00 S208 Oe 

TOTAL MEAN= 21.60 TOTAL S-DEV= 3.75 Y= 23 MAX-DMD= 66 

SIMULATION: AV E605 2—= 246.61 FOF SHORKtS= 0.9046 

SHSIMU= 30.36% SD-COST= Bees) | SD-#OF SH'tS= 0.0124 

NORMAL~APPROXIM.: AVE-COST= 238.36 #COF SHOK'S= 9.8938 

SUMMARY STATIS ores: 

SD/AVE= 1.36% PISIMA= 212.033 PINORA= 266.67$ 

DEVMON= -3.35% DEVNSH= -1.19% DEV?PLIA= -2.19% 

NOPFDMD= 3 #OF TIME UNITS= 72 REPLI=20000 I-SEED= 901941. 
#OF REP- PROBA- UNIT=-SHORT. rl 
CY CieEs Biola wy MEAN Cos N*P Hey 

1: 18 ie O20 le BOW lo 3100200 3.6. 1s 
ZO 2: 0.40 Ze 714.40 Ze 300.00 T24 2166 
3 2 le 327 0.30 32 5460 Be 100.90 3.6 Coe 

TOTAL MEAN= 21.60 TOTAL S-DEV= 3.75 Y= 20 MAX-DMD= 66 

SIMULATION: AV E-COST= 65.10 #OF SHOR! S= OezZ2 93 

SHSIMU= 9.80% SD-COST= 1.66 SD-#OF S#H'S= 0.0061 

NORMAL-~APPROXIM.: AVE-COST= 58.49 FOF SHOR'S= 0.2193 

SUUMARY SoA hoes. 

SD/AVE= 2-04% PISIMA= 213527295 PINORA= 266.67$ 

DEVMON= -7.36% DEVNSH= -4.37% DEVPIA= -3.13% 


APPENDIX D 


SIAULATION RESULTS FOR THREE DEMAND SOURCES. 


NOFDMD= 3 #OF TIME UNITS= 60 REPLI=20000 I-SEED= 1309 
40F REP-  PROBA- SNIT-SHORT. FIT 
CYCLES BILITY MEAN 705 Nx NXQ 

1: 10 1: 0.20 1; 2.00 i: 100.00 2.0 8.0 
2: 12 nO Dee 3.60 2: 300.00 3:6 «8.4 
3: 15 Senn ue 0 Be 12.00 3: 200.00 12.0 3.0 

SOTAL MEAN= 17.60 TOTAL S-DEV= 2.55 Y= 23  MAX-DND= ~37 

SIMULATION: AVE-COST= 49.11 #40F SHOR'S= 0.2333 

SHSIMU= 12.61% SD-COST= 1.09 SD-#0F SH'S= 010051 

NCRMAL-APPROXIM.: AVE-COST= 48.93 #0F SHOR'S= 0.2340 

SUMMARY STATISTICS: 

SD/AVE= 3.23% “PISIMA= 210.473 PINORA= 209.09$ 

DEVMON=  -0.37% DEVNSH= 0.29% DEVPTA= °-0.05% 

NOFDMD= 3 #OF TIME UNITS= 72 REPLI=20000 {I-SEED= 84009 
OF REP-  PROBA- UNIT-SHOP?. PIT 
CYCLES BILITY MEAN cos N*D  NXO 

mee Ai Ve Otho aes 4.80 1 300.00 Le ae 
2: 24 De EtG ; 442200 2: 300290 12.0 12:0 
Sc 3° 0.60 3: Or 60 3: 100.00 10:8 7.2 
TOTAL MEAN= 27.60 TOTAL S-DEV= 3. y= 32 MAX-DND= 54 
SIMULATION: AVE-COST= 45.55 #0OF SHOR'S= 0.2010 

SHSIMU= 9.05% SD-COST= 1.34 SD-#0F SH'S= 010054 
NORMAL-APPROXIM.: AVE-COST= 43.69 #0F SHOR'S= 0.1970 

SUMMARY STATISTICS: 

SD/AVE= 5.72% PISIMA= 226.57$ PINORA= 221.74$ 

DEVMON= -4.08% § DEVNSH= —-1.99% DEVPIA= -2. 13% 

NOFDMD= 3 #OF TIME UNITS= 60 REPLI=20000 I-SEED= 891. 
#OF REP- PROBA-_ UNIT-SHORT. FIT 
CYCLES BILITY MEAN Cost N+P N# 

de 1: 0.70 7: 10.50 1: 100.00 10.5 4. 
Di 20 2: 0230 ea O10 2: 300200 620 14.0 
3: 12 Seno Gio Bo S300 200 jp ee: 
TOTAL MEAN= 23.70 TOTAL S-DEV= 3.20 Y= 27 MAX-DMD= 47 
SIMULATION: AVE-COST= 49.95 40F SHOR'S= 0.2356 

SHSIMU= 11.55% SD-COStT= 1.32 SpD-#0F SH'S= 0.0055 
NORMAL-APPROXIM.: AVE-COST= 52.36 #OF SHOR'S= 0.2477 

SUMMARY STATISTICS: 

SD/AVE= 2.44% PISIMA= 211.973 PINORA= 211.39$ 

DEVMON= 4283%  ## DEVNSH= 5.11% DEVPIA= -0.27% 
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NOFDMD= 3 #OF TIME UNITS= 60 REPLI=20000 I-SEED= 312774 
#0F REP-  PROBA- UNIT-SHORT. PIT 
CYCLES BILITY MEAN cost N*D KO 

1: 10 j2m0ns 0 jeans 00 1: 300.00 BLO ero 
De 8 DO Ga >. Gao 2: 3 O0mon 370 6.0 
3: 20 SS 3: 10:00 3: 100.00 10.0 1020 

TOTAL MEAN= 22.00 TOTAL S-DEV= 3.27. Y= 27. £MAX-DMD= 45 

SIMULATION: AVE-COST= 17.29 #0F SHOR'S= 0.0904 

SHSIMUZ 4.88% SD-COST= 0:68 SD-#0F SH'S= 0.0033 

NORMAL-APPROXIM.: AVE-COST= 18.45 #0F SHOR'S= 0.0882 

SUMMARY STATISTICS: 

SD/AVE= -93%  PISIMA= 191.123 PINORA= 209.09$ 

DEV MON= "71%  ##DEVNSH=  —2.46% DEVPTA= 97.40% 

NOFDMD= 3 #OF TIME UNITS= 88 REPLI=20000 I-SEED= 932675. 
#0F REP-  PROBA- UNIT-SHORT. FIT 
CYCLES BILITY MEAN Cast *P xO 

1: 11 1: 9.10 1: 1.10 1: 100.09 1.1 9.9 
2: 44 2: 0:30 2: 13.20 2: 100:00 13:2 30.8 
aee? Sona 3 anc 00 3: 300.060 19 cee 

TOTAL MEAN= 34.10 TOTAL S-DEV= 3.49 Y= 39 MAX-DMD= 77 

SIMULATION: AVE-COST= 24.69 #9F SHOR'S= 0.1268 

SHSIMU= 5.86% SD-COST 0285 SD-#0F SH'S= 0.0043 

NORMAL-APPROXIM.: AVE-COST= 27.17 #0F SHOR'S= 0.1257 

SUMMARY STATISTICS: 

SD/AVE= 3.57% “PISIMA= 193.933 PINORA= 216.13$ 

DEVMON= 10.45% DEVNSH= -0.89% DEVPTA= 11.45% 

NOFDMD= 3 #OF TIME UNITS= 84 REPLI=20000 I-SEED= 928698 
40F REP- PROBA- UNIT-SHORT. FIT 
CYCLES BILITY MEAN cost N*P N*Q 

ao 1: 0.30 1: 6.30 1: 300.00 6. Sie 
2: 42 De OEO eee aod 2: 300299 25) 2 Es 
3: 14 2 nO See SO 3: 100200 523 °1112 
TOTAL MEAN= 34.30 ‘TOTAL S-DEV= 4.09 Y= 49 MAX-DMD= 77 
SIMULATION: AVE-COST= 39.55 #0F SHOR'S= 0.1366 

SHSIMU= 6.12% SD-COST= 1.31 SD-#0F SH'S= 0.9046 

NORMAL-APPROXIM.: AVE-COST= 42.72 #0OF SHOR'S= 0.1506 

SUMMARY STATISTICS: | 

SD/AVE= 3.32%  “pISIMA= 289.46$ PINORA= 283.67$ 

DEV MON= "00%  DEVNSH= 10.20% DEVPIA=  -2.00% 
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NOFDMD= 3 #OF TIME UNITS= 48 REPLI=20000 I-SZzD= 3311 
#OF REP- PROBA- UNI T- SHOR?. FIT 
CYCLES BILITY MEAN N*D  N#Q 
1: 12 1: 0.80 1: 9.60 4: 100.09 9.6 2.4 
2? 24 2: 0.40 2: 9260 2: 250290 $26 1424 
3: 16 33 0210 3: «8612360 «= 3: 250200 26 1424 
TOTAL MEAN= 20.80 TOTAL S-DEV= 3.02 y= 24 MAX-DMD= 52 
SIMULATION: AVE-COST= 43.50 #9F SHOR'S= 0.2196 
SHSIMUS 10-85% SD-COST= 1:05 SD-#0F SH'S= 0.0053 
NORMAL-APPROXIM.: AVE-COST= 39.98 #0F SHOR'S= 0.2212 
SUMMARY STATISTICS: 
SD/AVE= 2.42%  PISIMA= 198.10$ | PINORA= 180.773 
DEVMON= -8.09%  DEVNSH= © 0:72% ## DEVPIA= -8.75% 
NOFDMD= 3 #0F TIME UNITS= 80 REPLI=20000 I-SEED= 395113. 
#OF REP- PROBA- UNIT-SHORT. PIT 
CYCLES BILITY MEAN C Ne N*Q 
1: 20 1: 0.90 1: 18.00 1: 100.00 18.0 2.0 
2: 80 2.20 eset) Meee O00 «8©=— «24201 Son0 
3: 40 : 0:90 Se SeaT0 32 300800 3650 «al0 
TOTAL MEAN= 78.00 TOTAL S-DEV= 4.71 Y= 82 MAX-DMD= 140 
SIMULATION: AYE-COST= 102.12 #OF SHOR'S=_ 9.5370 
SHSIMU= 17.13% SD-COST= [98 SD-#0F SEH'S= 0.0105 
NORMAL-APPROXIM.: AVE-COST= 99.36 #0F SHOR'S= 0.5167 
SUMMARY STATISTICS: 
SD/AVE= 1.94% PISIMA= 190.178 | PINORA= 192.31$ 
DEVMON= -2.70% DEVNSH=  —-3:78% | DEVPIA= 1.12 
NOFDMD= 3 #OF TIME UNITS= 60 REPLI=20000 I-SEED= 7655599 
#0F REP- PROBA- UNIT-SHORT. IT 
CYCLES BILITY MEAN Cost N*D NKQ 
1: 10 1: 0.20 1: 2.00 14: 100.00 2.0 8.0 
2: 30 2: 0:40 Peuigato 2: | We00 122.0) 18n0 
32715 3: 0220 3: 3200 3: 300200 320 12:0 
TOTAL MEAN= 17.00 TOTAL S-DEV= 3.35 Y= 22 MAX-DMD= 55 
SIMULATION: AVE-COST= 14.71 #0F SHOR'S= 0.1114 
SHSIMU= 5.48% SD-COST= 0:55 SD-#0F SH'S= 0:0039 
NORMAL-APPROXIM.: AVE-COST= 13.23 #0F SHOR'S= 0.0973 
SUMMARY STATISTICS: 
SD/AVE= 3072% PISIMA= 131.97$  PINORA= 135.29$ 
DEVMON= -10.03%  DEVNSH= -12:24% DEVPIA= 2.52% 
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NOFDMD= 3 #OF TIME UNITS= 60 REPLI=2000) I-SEZD= 34876 
#9F REP- PROBA- UNI T-SHORT. FIT 
CYCLES BILITY MEAN cost *P NO 

1: 20 1: 0.20 1: 4.00 1: 300.00 4.0 16.0 
2: 30 2: 0:30 2: 9500 | 2: —BaGmo0 910 2120 
3: 75 3: 0:60 3: $100 = =§=©93: +=©100200 FO ace 

TOTAL MEAN= 22.00  OTAL S-DEV= 3.62 Y= 26 MAX-DUD= 65 

SIMULATION: AVE-COST= 57.84 #0F SHOR'S= 0.2477 

SHSIMU= 10.92% SD-COST= 1:44 SD-#0F SH*S= 0.0061 

NORMAL-APPROXIM.: AVE-COST= 53.10 #OF SHOR'S= 0.2434 

SUMMARY STATISTICS: 

SD/AVE= 20-49% PISIMA= 233.453 | PINORA= 218.133 

DEVMON= -8.20%  DEVNSH= “-1.77% | DEVPIA= -6.54% 

NOPDMD= 3  #OF TIME UNITS= 72 REPLI=20000 I-SEEZD= 176532 
#OF REP- PROBA- UNIT- SHORT. FIT 
CYCLES BILITY MEAN cosa *P N*Q 

1: 12 1: 0.90 1: 10.80 1: 100.00 10.8 1.2 
2: 36 2: 0:40 2: 18540 © 2: AI@ONOO 10.0 hime 
3: 718 3: 0:50 3: 9100 3: 300200 950 Seo 

TOTAL MEAN= 34.20 TOTAL S-DEV= 3.77 Y= 39 MAX-DMD= 66 

SIMULATION: AVE-COST= 26.99 #0F SHOR'S= 0.1781 

SHSIMU= 8.02% SD-COST= 0.80 SD-#0F SH*S= 0.0051 

NORMAL-APPROXIM.: AVE-COST= 27.45 #OF SHOR'S= 0.1799 

SUMMARY STATISTICS: 

SD/AVE= 2.96%  PISIMA= 151.55$  PINORA= 152. 633 

DEVMON= 1.71%  DEVNSH= 0.99%  #DEVPIA= —0:71% 

NOFDMD= 3 #OF TIME UNITS= 78 REPLI=20000 I-SEED= 3480610 
#0F REP- PROBA- UNIT-SHORT. FIT 
CYCLES BILITY MEAN COST N*D  N¥Q 

1: 13 1: 0.40 1: 5.20 1: 100.09 5.2 7.8 
2: 78 23 “ORO 2: 15160 2: igoeo0) Viscomemes 
3: 39 3: 0270 3: 27.30 32 SOUE008) 273mm 

TOTAL MEAN= 48.10 TOTAL S-DEV= 4.88 Y= 55 MAX-DMD= 130 

SIMULATION: AVE-COST= 37.89 #0F SHOR'S= 0.1828 

SHSIMU= 6.73% SD-COST= 1:23 SD-#0F SH*S= 010059 

NORMAL-APPROXIM.: AVE-COST= 36.71 #OF SHOR'S= 0.1719 

SUMMARY STATISTICS: 

SD/AVE= 3.25%  PISIMA= 207.273 | PINORA= 213.51 

DEVMON= -3.12%  DEVNSH= ~-5.96% DEVPIA= ~ 3.01% 


= ap ae com cee cee ee ee ee ce ee eo ee ee es es ee es es ee ee ee es ee ee =_— =P ab oe Se ee ee eee eee eee ee eee eee ee ee 
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NOFDND= 3 


#OF REP- PROBA- 
CYCLES Bags. ¥ MEAN 
We lt sO. 05 le 0.55 
2: 66 221102350 ve 19.80 
Bes eps, AO oa 3520 
TOTAL MEAN= 23.65 TOTAL S-DEV= 4 
SIMULATION: AVE-COST= 5.25 
SHSIMU= 5.46% SD-COST= 1.26 
NORMAL~APPROXIM.: AVE-COST= 31.08 
SUMMARY oLTATISTICS: 
SD/AVE= Be 5 Ti PISIMA= 272.73$ 
DEVMON= -11.34% DEVNSH= ~11.63% 


FOF TIME UNITS= 66 REPLI=20000 


NOFDMD= 3 #OF TIME UNITS= 72 REPLI=29000 


#OF REP- PROBA- 
CYCLES ee cele ny MEAN 
1 18 Tes O20 ae 3.60 
eer 36 2: 0.40 Ze 14.40 
ae. 12 32 0230 see Be oO 
TOTAL MEAN= 21.60 TOTAL S-DEV= 3. 
Si ATION: AVE-COST= 246.61 
SHS IMU= 30. ter SD-COS T= Siero, 
NORMAL~APPROXKIM. : AVE~COST= 2362510 


Puan s SLATESTICS; 
SD/AVE= 1.36% 
DEVMON= 


NOFDMD= 3 #OF TIME UNITS= 40 REPLI=20000 


tor REP= PRODA= 
CrCi ES BEE TY MEAN 
13 10 lee Uso As OU 
ZieeZ0 Zee a0) Cee 
4Q Bo U 5 ceed 
TOTAL MEAN= 35.00 TOTAL S-DEV= 4.11 
SUA T LON: AV E&=COST=° 122.37 
SHSIMU= 19.58% SUsCCS1— Zao 
NORMAL OAPPROXTM.  AVE-COST= 118.98 
Sane Sioa col LCS 3 
SD/AVE= Vet? PISIMA= 228.785 
BEVHON= -2.77% DsVNSH= Se 6 1% 


I-SEED= 394861 
UNIT-SHORT. 
N*D NEO 
1; 300.90 Onsen 
2: 300.00 19:8 46. 
3: 100.00 B2a52oen 
Y= 30 MAX=DMD= 119 
#OF SHOR'S= 0.1292 
SD-#0F SH'S= 0:9046 
#OF SHOR'S= 0.1142 
PINORA= 272.093 
DEVPIA= -0.23% 
I-SEED= 901941 
ONIT-SHORT. 
COST N*D  N*Q 
1: 300.00 seo 4 4 
2: 300.00 1454 2106 
3: 100.00 3.6) yee 
Y= 23 MAX-DMD= 66 
#OF SHOR'S= 0.9046 
SD-#0F SH'S= 0.0124 
#OF SHOR'S= 0.8938 
PINORA= 266.67$ 
DEVPIA= ~-2.19% 
I-SEED= 32467 
UNIT-SHORT. Ppa 
Soiam N*D ~NKQ 
1: 100.00 JO mare 0 
2: 100.00 1220 8.0 
3: 300.00 Djs ONO 
Y= 38  MAX-DMD= 79 
#0F SHOR'S= 9.5349 
SD-#OF SH'S= 0.0096 
#OF SHOR'S= 0.5552 
PINORA= 214.293 
DEVPTA= -6.34% 
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APPENDIX E 


FOR FIVE DEMAND SOUKCES 


—-— = qe com co co ce es se see es ee es es es es ee eee ee ooo oa eee ee eee eee eee ee ee eee ee ee ee 


REPLI=20000 


EV= 
333200 
a 


4.01 


SIMULATION RESULTS 
NOFDMD= 5 #OF TIME UNITS= 60 
TOP whe o> PROBA- 
CYCives ivi Ty MEAN 
1; 10 1; 0.10 i 
20 22 eue30 2° 
ae 3.2 50250 3: 1 
Geel Ye 0) 5748 Ye 
5 salto 5270290 5 1 
TOTAL MEAN= 41.90 TOTAL S-—D 
SI ATION: AV £E-COST= 
SASIMU= ioe nay SD-COST= 
NORMAL~APPROXIM.: AVE-COST= 
SUMMARY “SIATISTICS. 
SD/AVE= -09% PISIMA= Z 
DEVMON= 0.93% DEVNSH= 
NOFDMD= 5 #OF TIME UNITS= 60 
#OF REP- PROBA- 
CY Cis Bier Ty MEAN 
1: 10 le O96 1: 
2s 30 2 ea 0 oe 
ozo S50. 50 a 1 
GQ: 12 Ge 0.70 Y: 
5=. iS Feo 0 ae 1 
TOTAL MEAN= 41.90 TOTAL S=) 
Si ATION: AV E-COST= 
SHSIWNU= ie bats SD-COST= 
NORMAL-~APPROXIM. : AVE-COST= 


SUMMARY 
SD/AVE= 
DEVMON= 


2 aa cen 


2 0 


REPLI=20000 


E Y= 
69.44 
dz 


69.07 


4.01 


I-SEED= 243078 
UNITS sHORT: Eek 
COoe NSP NO 
ls) 20000 1.0) eed 
2: “300200 O's ON 2eln0 
3: TO0mG0 102.05 Foo 
4; 300.09 C24 woe 
33° 2UOROe 132:55 le 
Y= 42 TAX— DAD = ead), 
#OF SHOR'S= 1.5449 
SD-#0OF SHiis= 0.0 1G) 
FOF SHOR’ S— 1.5455 
PINORA= 217.66$ 
DEVPIA= 0. 89% 
J-SESD=)) 242076 
Get = SOK « Peat 
Soot N¥P N*¥Q 
lew 200 200 1. 0026 
2s) 300800 oO eno 
3: (100R20 POs Umeda) 
G4: 300.90 8. Gea 
9; 20000 list) Stes 
Y= 46 MAX-DMD= 8/7 
#OF SHOR'S= OlS225 
SD=#OF SH*S= 0.0074 
#OF SHOR'S= 02 3am 3 
PINORA= 217.663 
DEVPIA= 1.18% 


3Z 


NOFOMD= 5 #0F TINE UNITS="60  REPLI=2Z0000~~I=SEED="2ZA3973— 


#0F RED- PROBA- 
CYeCRES Beri Ty MEAN 
10 i Oe 10 1: OO 
230 Pea SU oe 9.00 
5 20 Bee Os 0) SF 10.00 
we 2 OR 718, ive Sree 
ite 1S be 10. 96 as ioe a0 
TOTAL MEAN= 41.90 TOTAL S-DEV= 4.01 
SIMULATION: AVE-COST= 24.09 
SHSIMU= 5.08% SD-COST= 0.92 
NORMAL-~APPECXIM.: AVE-COST= 2405 


SuunanY OLATISTICS; 
SD/AVE= Scola 
DEVMON= =Oehoy 


em ae em ee tt em mmm se se ee Se 
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UN PSOne. Boles 
COST Neaeaeeeh= O 
ie 200200 cee oO 
2: 300.00 Oe mee ler) 
oe 100.00 PO2 CeO 20 
4; 300.00 Sa 30 
S72 0200 13.5" aes 
Y= 48 4AX-DMD= 87 
#OF SHOR'S= Oiewiel SC 
SD-#OF SH'S= Q.0042 
#OF SHOR'S= Os (alos 
PINORA= 217.665 
DEVPIA= 221k 


—— ee a ce ee ee em i ee ee ee 


APPENDIX F 
OUTPUT TABLEAU OF THE SIMULATION PROGRAM FOR TWO DEMAND 
SOURCES 


He ke ke ok ee oe eo eK Re Ke KK ee 2 ee ee KK A OK KK KK KKK KK 


REPLIC*NS: 100000! TIME UNITS PER PERTOD:20-—) [-Siep i504 


#OPF REPCYCLS(1):4%  P1=0.255 MU1=1.00% PI1= 109.007 Y=78 
#OF REPCYCLS(2):59  P2=0.7519MU2=3.751!1PI2=1000.0012 


RESULT OF BERNOULLI TRIALS FOR BACAR REPAIR CYGiee 
ESTIMATED MEAN OF SHORTAGE COSTS PER PERIOD: 13.43733 
ESTIMATED ST=DEV. FORK ESTIMATES OF THE SNEAN 2 Sve oo 


RESULT OF MONOVAK. APPROXIMATING FORMULA: 11.771195 

RESULT OF HYPERGEO. APPROXIMATING FORMULA: 9.33216 

RESULT OF CDF APPROXIMATING FORMULA? Sie. 43 722 
SUMMARY:128MU= 4.7519 SHEXAC=1.36%20 SHSIMU= 1.36%21 


SD/MEAN=2.7%22 DEVMON=-12.4% DEVHYP=-30.6 % DEVCDF= 0.0%2 3 


MAE AR IK RK ee AK EK IK A A AK I FO KK I KK KKK EK OK 


1. Number of simulated periods; 

2- Length of one period; 

3. Initial seed, necessary to reproduce identical sinu- 
lation result; 

4. Number of repair cycles for end item no. 1 perc 
period. For instance, if the time period has a 


length of 20 time units and the number of fcepair 
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9. 
Vee 


1S 
lie 


ee 


14. 


lee 


V Oye 


Aye 


Tce 


132 


20. 


Zale 


cycles is 4 then each repair cycle has a length of 5 
time units. 

EGO rarity) SOpdchand Inweoemesrebaict cycle for end 
item no.1; 

dean demand for period for end item no.1; 

Unit shortage cost if one jemand caused tky end iten 
Now Not satisfied: 

Purchased quantity available to Satisfy period 
demand; 

Number of repair cycles for end item no.2; 
Probability of demand in one repair cycle for end 
item no.2; 

Mean demand for period for end item no.2; 

Unit shortage cost, one demand causei by end iten 
nho.2 not Satisfied; 

Sample mean of the shortage costs per period as 
observed during Simulation; 

Estimates precision of mean eStimate. 

Monovariable (McMasters') approximation formula; 
Hypergeometric approximation formula; 

CDF approximation fornula; 

The last two lines provide a condensed simulation 
Summary. 

Sum of expected demands from both demand sources for 
one period; 

The program computes the distribution of the sum of 
the two binomial demand distributions and is there- 
fore able to output an exact expected fraction of 
periods in which shortages occur. In this example 
1360 periods with shortages are the expected number 
so the expected fraction is 1.36%. 

Fraction of periods in which shortages occurred 
during Simulation. When compared to SHEXAC it 
provides a numerical hint for the quality of the 


random numbers generated during the simulation. 
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are 


23 


The ratio of estimated standard deviation for the 
estimate of the mean to the sample nmean. SD/MEAN 
decreases as the number of replications increases and 
the confidence interval for the mean derived from 
Simulation becomes smaller. The value of SD/MEAN is 
the basis for the decision as to whether the number 
of replications ensures sufficient precision. The 
number of replications must be traded off against the 
desired precision. 

DEVMON, DEVHYP, and DEVCDF express how £ar the three 
approximation formulas deviate from the Simulation 
resuit. The complete last row allows a guick survey 
about the performance of the approximations. Frog 
the observed sample mean for the shortage costs per 
period and its estimated standard deviation, a confi- 
dence interval for the true unknown mean can be 
computed. It can be used to determine the value of 
the confidence level such that the confidence 
inferval for the true mean would include the values 
of the approximation formulas. FOr LASstaticew a. 
SD/MEAN = 5% and DEVMON = -104, it follows that the 
value of the monovariable approximation formula is 
two estimated standard deviations away from the 
sample mean and would be covered by a confidence 
interval with a confidence level of at least .9545 
ise. P(IZ) <e2\e=s ooo 
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